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Abstract 

We study the ID extended Hubbard model with a weak repulsive short-range interaction 
in the non-half-filled band case, using non-perturbative Renormalization Group methods and 
Ward Identities. At the critical temperature, T — 0, the response functions have anomalous 
power-law decay with multiplicative logarithmic corrections. The critical exponents, the 
susceptibility and the Drude weight verify the universal Luttinger liquid relations. Borel 
summability and (a weak form of) Spin-Charge separation is established. 
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1 Main Results 



1.1 Introduction 

The Hubbard model, see e.g. [I], describing interacting spinning fermions on a lattice, plays the 
same role in quantum many body theory as the Ising model in classical statistical mechanics, 
that is it is the simplest model displaying many real world features: it is however much more 
difficult to analyze. While our understanding of the Hubbard model in higher dimensions at 
zero temperature is really poor (except for special choices of the lattice as in [2] and the 
situation is better in d = 1, when the model furnishes an accurate description of real systems, 
like quantum wires or carbon nanotubes [3]. 

The one dimensional Hubbard model (from now on the Hubbard model tout court) can 
be exactly solved by Bethe ansatz, as shown by Lieb and Wu [5]: the system is insulating in 
the half filled band case while it is a metal otherwise and the elementary excitations are not 
electronlike, a phenomenon which is nowadays called electron fractionalization [BJ. Recently in 
[7] a strategy for a proof that the lowest energy state of Bethe ansatz form is really the ground 
state has been outlined (see also [8]). This method is however of little utility for understanding 
the asymptotic behavior of correlations; and does not apply in studying the ground state of a 
slight generalization of the model, the extended Hubbard model, that consists in replacing the 
local quartic interaction with a short-ranged one. Other approaches has been therefore developed 
to get more insights into the physical properties of the Hubbard model. 

Under certain drastic approximations, like replacing the sinusoidal dispersion relation with a 
linear relativistic one and neglecting certain terms called backward and umklapp interactions (see 
after (|2.25[) for their definition) , one obtains the spinning Luttinger Model, which is exactly solv- 
able in a stronger sense, [5], |10) : all its Schwinger functions, at distinct points, can be explicitly 
computed. This model, describing interacting fermions, can be exactly mapped in a model of 
two kinds of free bosons, describing the propagation of charge or spin degrees of freedom and 
with different velocities (spin-charge separation); again, a phenomenon of electron fractional- 
ization which has received a considerable attention in the context of high T c superconductors 
Moreover, as in spinless Luttinger model, the correlations have a power law decay rate with 
interaction dependent exponents. 

However, neglecting the lattice effects and backscattering or umklapp interactions is not 
safe, and indeed the mapping to free bosons is not expected to be true in the Hubbard model. 
A somewhat more realistic effective description can be obtained by including the backward 
interaction in the spinning Luttinger Model, so obtaining the g-ology model. This system is no 
more solvable; however, a perturbative Renormalization Group (RG) analysis, |12) . shows that, 
in the repulsive case, such extra coupling is marginally irrelevant, i.e. becomes negligible over 
large space-time scales. In [13] the necessity of implementing Ward Identities in a RG approach 
was emphasized, in order to go beyond purely perturbative results, but the analysis was limited 
to the Luttinger model and no attempt was done to include the effects of nonlinear bands. In 
[14] it was observed that the correlations of the repulsive g-ology model would qualitatively differ 
from the Luttinger model ones for the presence of multiplicative logarithmic corrections. 

A new point of view, that extended previous ideas of Kadanoff, [15) . and Luther and Peschel, 
[16j . was proposed by Haldane, [17] . and is nowadays known as Luttinger Liquid Conjecture. The 
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idea is to exploit the concept of universality, a basic notion in statistical physics saying that the 
critical properties are largely independent from the details of the model, at least inside a certain 
class of models. In the present the exponents are non trivial functions of the coupling, 

universality has a meaning more subtle than usual; it does not mean that the exponents are the 
same (the exponents do depend on the details of the model), but that the exponents and certain 
thermodynamic quantities verify a set of universal relations which are identical to the Luttinger 
model ones. Such relations give an exact determination of physical quantities in terms of a 
few measurable parameters. The validity of such relations has been checked in certain special 
solvable spinless fermionic lattice models |17j , but a proof of their validity in the Hubbard model 
(or in the non solvable extended Hubbbard model) is an open problem. It should be remarked 
that, even thought the Hubbard model differs from the spinning Luttinger model for irrelevant 
or marginally irrelevant terms in the Renormalization Group sense (in the weak non half filled 
band case), this would not imply at all the validity of the same relations as in the Luttinger 
model; irrelevant terms do renormalize the exponents and the thermodynamic quantities. 

Starting from the 90's, the methods developed in constructive Quantum Field Theory (QFT) 
for the analysis of QFT models in d = 1 + 1 [HI [H] were applied to interacting non relativistic 
spinless fermionic systems in the continuum |20) . so establishing the anomalous dimension in a 
non solvable model, by combining Renormalization Group methods with non perturbative infor- 
mation coming from the exact solution of the Luttinger model] the extension to spinning fermions 
was done in [21j . An important technical advance was achieved in [221123] , by implementing Ward 
Identities based on local symmetries with Renormalization Group methods. A well known diffi- 
culty in any Wilsonian Rermalization Group approach is that the momentum cut-offs break the 
local symmetries on which Ward Identities are based; in [22l [23] it was developed a technique 
allowing to rigorously take into account the extra terms produced in the Ward Identities by the 
cut-offs, so that interacting non relativistic fermions in d = 1 were constructed without any use 
of exact solutions [231 1241 [2"5] . The main outcome, with respect to the early Renormalization 
Group analysis [J2] , is that the results are exact (the lattice and non linear bands are fully taking 
into account) and non-perturbative; and physical observables are written in terms of convergent 
expansions so that they can be computed with arbitrary precision. The complexity of such ex- 
pansions made however impossible to verify explicitly the universal Luttinger Liquid relations 
in J2S]; they have finally been proven for interacting spinless fermions on the lattice (the XX Z 
chain and extended versions) in [2H] through Ward identities; this fact appears to be related 
to a non-perturbative version of the Adler-Bardeen theorem of the non renormalization of the 
anomalies, .29 . In this paper we will extend such ideas to spinning fermions in the Hubbard 
model; as we will see, the extension is rather non trivial and new phenomena take place. 

1.2 Extended Hubbard Model and Physical Observables 

Let /3 > be the inverse temperature, —[/, the chemical potential and C — {—[L/2], . . . , |[(L — 
l)/2]L} a one dimensional lattice of L sites. The extended Hubbard model [30 describes fermions 
hopping on C with a short-range density-density interaction; the Hamiltonian plus the chemical 
potential term is 

H = -\/Z( a ts a x+l,s + a t,s a x-l,s) + V/Z a ts a x,s + X Yj V ( X -y) a t,s a x,s a ts> a y,s' ( L1 ) 

x£C x€C x.VQC 

s = ± s = ±l s,s'=±l 

where s are fermionic creation and annihilation operators at site x with 'spin' s, and v(x) is 
a smooth, even potential such that |i>(x)| < Ce _K ' x ' (short range condition); periodic boundary 
conditions are assumed: a £ L+1 s = a\ s . 

The set-up of the Grand Canonical Ensemble is standard, and we remind it concisely; more 
details are, for example, in [31j . If O x is a monomial in the operators a e x s , e,s — ±, or in 
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the density operators a% s a% s ,, s,s',s,s' = ±, given xq e [0,/3], define x = (x,xq) and O x := 
e Hx °O x e~ Hx ° (so that xo has the meaning of imaginary time); then, given the observables 
Xl , . . . , O x „ , their Grand Canonical correlation is 



Tr[e- ^T(Q X1 ■ ■ ■ Q x J] 
Tr[e-P H ] 



(O -O ) — — ^- 1 1 " 7J (12) 



where T is the time order product. Similarly, (0 Xl ;--- ',O x „) T . L g denotes the corresponding 
truncated correlations. We are interested in the correlations in the thermodynamic limit L — > oo 
and at the critical temperature = 0; the limit L,j3 — > oo will be indicated by dropping the 
labels L, f3. 

Define pp G [0,7r], the free Fermi momentum, and Vp, the free Fermi velocity, such that 

cospF = fi vp = sinpp . 

In this paper we have three main assumptions on the parameters: 

pp ^ 0, vr/2, tt , v(2p F ) > , A > (1.3) 

The condition pp ^ 0,7r means that the empty band and the filled band cases are not included; 
the reason of such exclusion is that, if vp — 0, the scaling of the model would be very different 
and would depend in a critical way on the interaction. The condition pp ^ it/2 excludes the 
half-filled band case; it will have the effect to make the Umklapp interaction (see fl2.ip irrelevant 
(in the RG language) . The two other conditions can be loosely called the repulsive condition on 
the interaction; they indeed imply that one of the contribution to the effective interaction (in 
the RG language) is strictly positive at all scales. 

The model is SU(2) symmetric, as the Hamiltonian is invariant under transformation a^ s — > 
J2 S ' Ms,s'0^ s > with M £ SU(2); and includes the standard and the U-V Hubbard models, 
corresponding to the interactions Xv(x — y) = US x . y and Xv(x — y) — U5 XtV + ^VS\ x ^ y \,i, 
respectively: in the former case the repulsive condition is U > 0. 

By definition, the T — free energy is 

E(X):=- lim (L^logTrte-^] , (1.4) 

and the 2-point Schwinger function is 

S-z,p,l{x - y) := (a~ + a+ + )p tL = (a~_a+_) 0<L . (1.5) 

The connection with experimental physics is through the response functions, defined as Fourier 
transforms of the following truncated correlations: 

JW(x - y) == (p>y) T ^, L ■= (p>y)(,, L - (f£)p,M}fi,L C 1 - 6 ) 

where p x is one of the following densities (see pagg. 54, 55 of [1]): 



a x . s a x,s (charge density) 

s=± 

Px = a x,s a( S a *,s' (spin densities) 

8,B'=± (1.7) 



- sa x,a a x,-a (singlet Cooper density) 

e = ± 

— ^ °x,«^»*«'°x+e,a' > e= (l)0) (triplet Cooper densities) 



»'=± 
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where i = 1,2,3 and 



" =' v o oj ^ = yi o y a =yo i 

When the interaction is off (A = 0), all functions 12 Q (x — y) have power law decay for large |x — y| 
with the same exponent; as we will see, turning on the interaction removes such degeneracy: some 
correlations are enhanced by the presence of interaction and others are depressed, so that the 
exponents become different. 

Define the Fourier transform of the response functions as 



n a (p):= lim / dioV^^Lfx) (1.8) 

where p = (p,Pa), with p e and po G The susceptibility is given by0 

k := limOcO.O) . (1.9) 

p— ¥0 

The paramagnetic part of the current J x is defined as 

^ x = o7 / y [°z+l.a a g.» _ a z,s a :r + l,s] (1T0) 
s=± 

while the diamagnetic part is 

Tx = ~\ \. a t,s a x+l,s + a x+l,s a x,s\ (1-11) 

The Drude weight is defined as 

D = -(t x ) - lim hm lim / dx } e lpx (J x J ) T . T B = lim lim Dtp) (1.12) 

po-S-0 p->-0 /3,Z,-Kx> ^ A ' po^Op^O 

where the first term is a constant independent from x. If one assumes analytic continuation in 
Po around po = 0, one can compute the conductivity in the linear response approximation by 

the Kubo formula, that is a = lim^^o h m 5->-o Z^+j ■ Therefore, a nonvanishing D indicates 
infinite conductivity. 
The conservation law 

^ = e H *°[H,p x ]e- Hx ° = -idPj x = -i[J x , Xa - J x -i, X0 ] , (1.13) 
dx 

where d x ^ denotes the lattice derivative, implies exact relations, called Ward identities (WI), 
between the Schwinger functions, the density correlations and the vertex functions, defined as 
Gp'-^x, y, z) = (pi C 'a y a+) T and G^ 1 (x, y,z) = (J x eiya+) T . Some Ward Identities, which will 
play an important role in the following, are 

-ipoG^Qc, k + p) - i(l - e-^Gf'ik, k + p) = S 2 (k) - 5 2 (k + p) (1.14) 

-z Po r!c(p) - t(l - e- i f)0 J> (p) = (1.15) 

-ip £l Pjj (p) - i(l - e~* )£»(p) = (1.16) 

where O i)(9 (x,y) = (p%Jy) Tt/3tL - 



1 in a fermion system, ft = K c p 2 , where k c is the compressibility and p the density, see e.g. (2.83) of 1321 
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1.3 Anomalous exponents and logarithmic corrections 



In the free case, A = 0, pf determines the position of the two singularities of the Fourier transform 
of 5*2, which are at k = (±j>f,0); whereas vf is the velocity of the large-distance leading term 
ofS 2 , 



S 2 (x) ~ Yl 



u=± 



In the following theorem we see that, when the interaction is turned on, A > 0, the singularities 
of the Fourier transform of £2 are moved into new positions, k = (±pp, 0), with pp = Pf + O(X). 
Moreover, the power law decay of S^x) and the response functions is strongly modified: the 
decay exponent is changed (anomalous dimension) and there are logarithmic corrections. 

Theorem 1.1 If the conditions \1.3\) are satisfied, there exist Ao > such that, i/0 < A < Ao, 
there exist continuous functions 

p F = pf(A) = Pf + 0(A) v F = u.f(A) = s'mp F (X) + 0(A) 

(depending also of the other parameters of the model, like v(x) and \x), such that, setting 



x := (x, vpxo) , L(x) = 1 + b\v(2pp) log |x| , b — 2(n sinp^)" 



n (x) 



5 (x) 



VpXQ COSPf — XSllipF 



(1.17) 



the large |x| asymptotics of the two-points Schwinger function is 

S 2 (x)~ [S (x) + i? 2 (x)] i(Xi " 



(1.18) 



where i?2(x) is a continuous function of X andx, such that |i?2(x)| < C^A 1 , for some positive 
constants 0# and 1? < 1. Besides, the large |x| asymptotics of the correlations are 



for a — C,Si 
for a = SC 
for a = TC, 



fi„(x) 
fi Q (x) 
Oq,(x) 



O (x) + i? a (x) , 1 £(x)C 
+ cos 2p F z 

L(x)£> 



fio(x) + i? Q (x) cos(2p F :r) 



l + i? Q (x) j 
1 L(x)t 



L(x)<° 



[l + i?a(x)] 



[l + i? Q (x)] 
(1.19) 



im'i/l £/ie functions R a (x) and i? Q (x) having the same properties of R2 (x) . Moreover, the critical 
exponents rj and X a , are continuous functions of X, while the exponents (sc and Ca of the 
logarithmic corrections could also depend on x (we can not exclude it), but satisfy the bounds 
\Csc\ < OA and |£ Q — £ Q | < CA, for a suitable constant C, with 



c. = , <c 



1 



(1.20) 



In the free A = case the response functions decay for large distance with power laws of 
exponent equal to 2. The interaction partially removes such degeneracy by producing anomalous 
exponents which are (in general) non trivial functions of the coupling (see Theorem 1.2); in 
particular the response to charge and spin densities are enhanced by the interaction, while the 
response to triplet Cooper densities are depressed. While the presence of non universal exponents 
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is a common feature with the Luttinger model, the presence of logarithmic corrections is a striking 
difference. Such corrections remove the degeneracy in the response of charge and spin densities: 
the response to spin density is dominating. Note on the other hand that the exponents of the 
non oscillating part of charge or spin density correlations are the same as in the free case; also 
logarithmic corrections are excluded. 

Similar formulas have been derived in the physical literature |14| under the g-ology approxima- 
tion, that is replacing the Hubbard model with the continuum g-ology model describing fermions 
with linear dispersion relation. The existence of anomalous exponents was proved previously in 
Theorem 1 of (51] (see also |H]), and the absence of logarithmic corrections in the non oscillating 
part of charge or spin density correlations was also previously proved in [33J . The above theorem 
improves such results, as it proves for the first time the existence of logarithmic corrections and 
universal relations. 

1.4 The Luttinger liquid relations 

Theorem 1.2 Under the same conditions of Theorem ll.ll there exist continuous functions 

K = K(X) = l-cX + 0(X 3/2 ), K = K(X) = l-cX + 0(X 3/2 ) (1.21) 

with c = [£(0) — v(2pf)](it sinpp ) , such that the critical exponents satisfy the extended scaling 
formulas 

Ari = K + K~ 1 -2, 2X C = 2X S . = K + 1 , 

(1.22) 

2X TCt = 2X SC =K~ 1 + 1, 2X SC = K + K . 



Moreover 



K v 2 p 2 
irv Pq + v A p A 

D(p) = -K P\ 2 +B(p) 

IT Pq + V p 



(1.23) 



with A(p), B(jp) continuous and vanishing at p = 0, v — sinpp + O(X) and K = 1 + 0(A); 
therefore the Drude weight D ll.lefjl and the susceptibility n hl.9\) are 0(A) close to their free 
values and verify the Luttinger liquid relation 

v 2 = D/k (1.24) 

The above Theorem says that, even if the logarithmic corrections alter the power law decay 
of the spinning Luttinger model, the exponents verify the same universal relations (ll.22[) . in 
agreement with the Luttinger liquid conjecture |15l 1161 117] . Such relations say that the knowledge 
of a single exponents implies the determination of all the others. 

The Fourier transform of the density correlation is similar to the free one, the interaction 
producing a renormalization of the velocity and of the amplitude K. The susceptibility and 
the Drude weight are finite, saying that the system has a metallic behavior (contrary to what 
happens in the half- filled band case). 

Besides the universal relations involving the critical exponents, there is also the universal 
relation (|1.24p . which relates the susceptibility and the Drude weight to the charge velocity v 
appearing in f| 1 . 23|) : this relation was conjectured in [T7] (in the spinless case, but the extension 
to the spinning case is straightforward, see e.g. g]). In the notation of [T7], vn — itk^ 1 ,v.j = ^ 
so that (]1.24j) takes the form 

v N vj = v 2 (1.25) 
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The validity of (jl.22p . (jl.24D is a rather remarkable universality property following from the 
combination of conservation laws of the Hubbard model and Ward Identities coming from the 
asymptotic gauge invariance of the effective theory. Wether a similar relation holds also for the 
spin conductivity is an interesting open problem. 



1.5 Spin Charge separation 

Theorem 1.3 Under the same condition of Theorem ll.li the Fourier transform of the 2-point 
Schwinger function is given by 

£ 2 (k + p£) = £(k)S M)W (k)[l + ii!(k)] , p F = (up F ,0) (1.26) 

where 

Z(k) = L(\k\- 1 )^[l + R'(k)}, |#(k)| < C|A| (1.28) 

L(t), t > 1, is the function defined in Theorem and Sm,u(M) * s a function whose Fourier 
transform is of the form 

Sm, u x = [p 1 F) J . ^^+o(i/|x|) L29 

ZtTVf \VpXo + IUJXi/Vf) ' [VcrXQ + IUJX\ / Vp ) 1 

with Vp t(T = 1 + 0(A), rjp = 0(A 2 ), v p — v a — c v \ + 0(X 2 ), with c v ^ 0. 

The above theorem says that the two point function can be written, up to a logarithmic 
correction, as the 2-point function of the spinning Luttinger model, a model which shows the 
phenomenon of spin-charge separation (see also |25|). A manifestation of spin charge separation 
is that the 2-point function is factorized in the product of two functions, similar to Schwinger 
functions of particles with different velocities. In this sense, the above theorem says that the 
spin-charge separation occurs approximately also in the Hubbard model, but is valid only at 
large distances and up to logarithmic corrections. Similar expressions are true also for the 
density correlations (the explicit formulae are in §5 and are not reported here for brevity). In 
the spinning Luttinger model v p = v, where v is the velocity appearing in (jl.23p ; in the present 
case we can verify this identity only at the lowest order in A, and whether this identity holds or 
not in the Hubbard model is an interesting open problem. 

1.6 Borel summability 

In ij2.6l we shall prove the following Theorem. 

Theorem 1.4 Given 5 £ (0,7r/2), there exists e = e(S) > 0, such that the free energy, the 
Schwinger functions and the density correlations are analytic in the set 

D e<5 = {A e C : |A| < e, \Arg A| < tt - 6} (1.30) 

and continuous in the closure, D £ ^$. Moreover, *//(A) is one of these functions, there exist three 
constants cq, c\, k, and a family of functions //i(A), h < 0, analytic in the set 

Co 



(1.31) 



\h\. 
such that 

o 

|A(A)| < de-^l , /(A)= £ f h (X) (1.32) 



h=—oo 



By using the Lemma in 119] . see pag. 466, this Theorem implies that all the functions satisfy 
the Watson Theorem, see pag. 192 of (34|. Hence they are Borel summable in the usual meaning. 
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1.7 Contents of the paper 

The paper is organized in the following way. 

1 . In i i2.1H2.3l we resume the RG analysis of the extended Hubbard model given in [35j [24] • 
The fermionic field is decomposed as a sum of fields ip( h \ h integer and < 1. The field 
ip^ 1 ' is associated to the momenta far from the Fermi points, while the fields ipw with 
h < are supported closer and closer to the two Fermi points, hence are more and more 
singular in the infrared region. The iterative integration of such fields, accompanied by 
a free measure renormalization (the field strength renormalization) , leads to a sequence 
of effective potentials V^ h \ expressed as renormalized expansions in a set of 5 running 
couplings Vh = (yh, Sh, 9i,h, g2,h, 9&,h), whose flow is driven by a recursive relation called 
beta function. The coupling Vh is associated to the only relevant term (in the usual RG 
language) present in the effective potentials (the others are marginal) and describes the 
change of the Fermi momentum due to the interaction; in order to control its flow, we 
change the chemical potential \x in /i — v and we compensate this operation by adding to 
the interaction a counterterm v s a+ s a~ s . The value of v is then chosen, by an iterative 
argument, so that Vh as h —> — oo; this will implicitly determine the interacting Fermi 
momentum. This procedure works since one can prove, under the conditions of Theorem 
11.11 the convergence of expansions in the running couplings, by using two crucial technical 
tools: the determinant bounds for fermionic truncated expectations and the partial vanish- 
ing of the beta function (see (I2.45[) below for the definition), which implies the convergence 
of Vh to V-oc = (0, S-oo, 0, 02,-00) 04,oo )< This limit can be seen as characterizing a point in 
a set of fixed points, depending on 3 parameters, of the RG transformation, suitably scaled; 
of course, the chosen fixed point depends on all details of the model. 

2. In Appendix [X] a property of the effective couplings flow is proved, implying the conver- 
gence of our expansions in the set D e< $, defined in (|1.30[) . At this point, as explained in 
£ 12.61 simple dimensional arguments allow us to prove Theorem 11.41 and then the Borel 
summability of the free energy and all the correlation functions. 

3. In £)2.4H2.5l the analysis is extended to the 2-point function and to the response functions, 
which are expressed in terms also of renormalizations of the density operators. In partic- 
ular, from the flow of these renormalizations one obtains the critical exponents and the 
logarithmic corrections appearing in (|1.19p . The critical exponents only depend oiidj? and 
v-oo, what will play a crucial role in the subsequent analysis. This is apparently not true 
for the logarithmic corrections; such corrections, absent in the spinless case, are due to the 
weaker convergence of the effective couplings to their limiting values. 

4. The Luttinger liquid relations (| 1 . 22[) or (|1.23[) can be checked directly by the expansions 
at lowest orders, but the complexity of such expansions makes essentially impossible their 
proof at any order. Similarly, the partial vanishing of the beta function, on which the 
RG analysis is based, cannot be proven directly from the expansions. Such properties are 
related to the asymptotic validity of certain symmetries, and our strategy consists in the 
introduction of a suitable effective model, for which such symmetries are exact, and in 
showing that certain quantities computed in the effective model coincides, with a proper 
choice of its parameters, with analogues quantities in the extended Hubbard model. The 
effective model is introduced in £!3.1l and is expressed directly in terms of functional integrals 
with linear dispersion relation; the model has an infrared and an ultraviolet momentum 
cut-off, and several interactions are present, non local and short ranged (both in space 
and time). The model can be consider a variation of the g-ology models introduced in 
the physical literature, the main difference being that the cut-offs are on space and time 
momentum components, what is of advantage for our approach (but makes the model not 
accessible to bosonization techniques). 
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5. The non locality of the interaction allows us the removal of the ultraviolet cut-off and a 
Renormalizaton Group analysis in the infrared region can be performed (similar to that 
performed in the Hubbard model), leading to convergent expansions in the effective cou- 
plings, see Appendix [B] A first use of the effective model is in the proof of the partial 
vanishing of the beta function of the extended Hubbard model; a proof of this property 
was already given in [24] . but we present here a simplified version of it (the main novelty 
is that the ultraviolet cut-off in the effective model is removed) , see Appendix [Cj Ward 
Identities and Schwinger-Dyson equations, with corrections due to the infrared cut-off, 
can be combined in the absence of back-scattering interactions to get relations implying 
the partial vanishing of the Beta function of the effective model; from this fact and using 
the symmetry properties in Appendix [Bj we can derive the partial vanishing of the beta 
function of the extended Hubbard model. Note the remarkable fact that a model with 
no back-scattering interaction and not spin-symmetric is used to prove properties of the 
Hubbard model, which is spin symmetric and in which the back-scattering interaction is 
present. 

6. If the back-scattering coupling is set equal to zero and both the infrared and ultraviolet cut- 
offs are removed, the model becomes exactly solvable, in the sense that the Schwinger func- 
tions verify a set of closed equations, obtained combining Ward Identities and Schwinger- 
Dyson equations, derived in !i3.2f43.5L The non locality of the interaction has the effect 
that the anomalies in the Ward Identities verify the Adler-Bardeen non renormalization 
property, see |36| . [28] . so that they can be exactly computed; therefore also the critical 
exponents, which are expressed in terms of such anomalies, can be exactly computed and 
the analogue of the relations (|1.22D for the effective model are obtained. 

7. In ^3.61 we prove that the limiting values V-oo of the effective couplings of the effective 
model with no back-scattering coincide with that of the Hubbard model, if a suitable fine 
tuning of the bare couplings of the effective model is done. Since the critical exponents 
only depend on V-oo and have the same functional dependence on it as in the Hubbard 
model, this implies that also the critical exponents coincide. Therefore we can prove that 
(|1.19p is satisfied for the extended Hubbard model, with the critical exponents verifying 
the relations (Tj~2"!Z|) . 

8. In <g]the proof of Theorem 11.31 is presented. By using the closed equation of the effective 
model in the limit of removed cut-offs, we can prove for it the exact Spin-Charge separa- 
tion. We show that this implies the approximate Spin-Charge separation for the extended 
Hubbard model. 

9. Finally, in $5] we compute the Drude weight and the susceptibility and we prove the Lut- 
tinger liquid relation (|1.24l) . The proof is based on the fact that these quantities are related 
to the Fourier transform of some density correlations. However, the bounds obtained in 
Theorem 11.11 in coordinate space do not allow us to exclude logarithmic singularities. In 
order to prove their finiteness, we use the Ward Identities of the Hubbard model (|1.14p . 
(ll.lSp . (|1.16[) . combined with the information coming from the effective model, keeping in 
this case also the backward interactions; this works since, even if the effective model is not 
completely solvable in that case, the Fourier transforms of the density correlations can be 
still exactly computed from the Ward Identities. We can prove that, by a suitable fine 
tuning of the couplings of the effective model, the Fourier transforms of the current and 
density operators coincide up to a constant and a renormalization, whose values are fixed 
by the Ward Identities (fTTT4"l) . ([TTT5]) . (TTTB)) : this implies (TOI)) . 



10 



2 RG Analysis for the Hubbard Model 



2.1 Functional integral representation 

The analysis of the Hubbard model correlations is done by a rigorous implementation of the RG 
techniques. To begin with, we need a functional integral representation of the model, because 
the RG techniques are optimized for that. We give here a concise description of it; a thorough 
discussion is in Sec 2 of [23]. The main object to study is the functional W( J, rj) = Wm-,l,p(J, v)i 
defined by 

e W(J,77) = J p^y-vm+J2 a fd^pZ+j: s fd x [ v l^- s +^+^- a ] ^ ^ 

where ipx s ano - Vx s are Grassmann variables and the fermionic density operators are de- 
fined as in (jl.7p . with i/j^s m place of s , J£ are commuting variables, J dx is a short form 

f° r SxgC /— ,a/2 d%o, P(dip) is a Grassmann Gaussian measure in the field variables tp^ s with 
covariance (the free propagator) given by 

P(#) iP^yy = , 1 P(#) i/>-,il>+- a = , 

In the above formulae, x(t) is a smooth compact support function equal to 1 for \t\ < 1 and 
equal to if \t\ > 7, for a given scaling parameter 7 > 1, fixed throughout the paper; M is a 
positive integer; £> i)/9 := V L x 2?^, £> L := %C, 2?,g := 2*(Z + \)\ 

vty) = A E / dxd y <^x,.«(x - y^+Ay,* ( 2 - 3 ) 

with w(x — y) = <5(a; — yo)v(x — y). Due to the presence of the ultraviolet cut-off r ) M , the 
Grassmann integral has a finite number of degree of freedom, hence it is well defined. The time 
shift in (12.21) . 8m '■= j3/y~M, is introduced in order to take correctly into account the discontinuity 
of g(x) at x = 0: our definition guarantees that, fixed L and /?, liiriM->oo <?m(x) = .9( x ) for x ^ 0, 
while limAf^oo <7a/(0, 0) = g(0, 0~), as it is to be for Proposition 12.11 below. 

If A = 0, the Hubbard model correlations can be easily calculated by using (|2.2[) . hence they 
are singular at momenta (cjpj?,0), uj = ±1. Since in the interacting theory, A 7^ 0, the position 
of the singularity is expected to change by O(A), when the first of conditions (|1.3j) is satisfied, 
we add to the interaction a counterterm 



s=± J 



and, to leave unchanged W( J, n) in (I2.1[) . we subtract the same term from the free measure, that 
has then a covariance: 

J_ v X ( 7 - M k )e^e~^ 
9M[ > ~ 0L ^ -ik + (cos PF - cosfc) ' 1 ' 

where pf is the interacting Fermi momentum defined such that 

cos pf = M — v . 
We introduce the following Grassmann integrals: 

5f^ L (x 1 , 8l , ei ;....;x n , Sn ,e n ) = - — 9 " _ £ W(J,n) (2.5) 
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It is well known that such Grassmann integrals, called Schwinger functions, can be used to 
compute the thermodynamical properties of the model with Hamiltonian (11.11) . This follows 
from the following proposition. 

Proposition 2.1 For any finite (3 and L, there exists a complex disc, centered in the origin, 
DL t p, such that, if X £ -Dl,/3, 

rr r -am ; ^Ht = Jim S n ,p ' (xi,si,£ 1 ;....;x n ,s n ,£ n ) ; (2.6) 

lr|e 1 \ M— ¥oo 

besides both members are analytic in X in the same disc. A similar statement holds for the 
density correlations. 

The proof of this theorem can be done exactly as in the spinless case, see |35j . The main 
point, strictly related with the fact that we are treating a fermionic problem, is that, for L 
and M finite, the l.h.s. of (|2.6p is the ratio of the traces of two matrices whose coefficients 
are entire functions of A, hence it is the ratio of two entire functions of A. Then, it may have 
a singularity only if Tr[e~^ H ] vanishes, which certainly does not happen in a neighborhood of 
A = small enough. On the other hand, it is rather easy to prove that also the r.h.s. of (|2.6[) 
is analytic in a small neighborhood of A = and that its Taylor coefficients coincide with those 
of the l.h.s.. This follows from the fact that the UV singularity of the free propagator is very 
mild and can be controlled with a trivial resummation, in the RG expansion, of the tadpole 
terms (see pag. 1383 of [35]). In this resummation, the only important thing to check is that 
limM-i-oo 3m(0, 0) = g(Q, 0~) (otherwise the perturbative expansions of the two sides of f|2 . 6[) 
would not coincide). 

The RG analysis will allow to prove that the analyticity domain is indeed of the form Dp ^ = 
{A,|A| < ceominKlog^rMlogL)- 1 }}^! < £ ,|argA| < § +5}, with c,e > 0, < 6 <ir/2 
independent of (3 and L. 



2.2 Multiscale analysis for the effective potential 

We will briefly recall here the RG analysis for interacting fermionic systems on the lattice as 
developed in [35] and [23] in the spinless and the spinning case, respectively. Note that the 
proof of many technical points do not depend on the spin, hence we shall refer to |35j for the 
corresponding details. 

Let T be the one dimensional torus, ||fc — k'\\j the usual distance between k and k' in T and 
||fc|| T = \\k - 0|| T - We introduce a positive function x(k') € C°°(T x R), k' = (k Q ,k'), such 
that x(k') = x( — k') = 1 if |k'| < to = aovp/^ and = if |k'| > a^vp, where vp = smpp, 
a = min{^, ™~% F } and |k'| = \Jk^+ v F \\k'W^. The above definition is such that the supports 
of x(k — Pf, ko) and x(fc + pf, ko) are disjoint and the C°° function onTxiJ 

A(k) := 1 - X (k - PF, k Q ) - x(k+PF, M (2.7) 
is equal to 0, if v F \\ [\k\ — pp] ||x + k^ < ig. We define also, for any integer h < 0, 

A(k')-x(7^k')-x(7- ,l+1 k') (2.8) 
which has support toj < |k'| < io7 /l+1 an d equals 1 at |k'| = toj , then 

X(k')= J2 M k ') ( 2 - 9 ) 

h=h,L,p 

where 

h L ,p := min {h : i o7 ' l+1 > |k m |} for k m = (tt/P,tt/L) . (2.10) 
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For h < we also define 

A(k) = f h (k-p F ,k ) + f h (k+p F ,k ) (2.11) 

(for h — 1 the definition is (|2.7[) ). This definition implies that, if h < 0, the support of /^(k) is the 
union of two disjoint sets, and A~^. In A^, k is strictly positive and ||fe — pf\\t < hi' 1 < io, 
while, in Aj^, k is strictly negative and \\k + p F \\j < ^ot' 1 - The label h is called the scaie or 
frequency label. Note that 



1= ]T A(k) (2.12) 



h—hr 



and that, since p F is not uniquely defined at finite volume (we are interested only to the zero 
temperature limit), then we can redefine it as (2ir/ L)(n F + 1/2), with n F — [Lp F /(2tt)]. Hence, 
if V' L = ?f(C + |) and V' L f} =V' L xV /3 ,we can write: 

.g(x-y)=.g( 1 )(x-y)+^ ]T e^'^-v)^)^ - y) (2.13) 



i,/3 



where 



9 [ y} ~ PL ^ e -iko + lcospp-cosk) [ ' 

5 £' i) (x-y) = ^ r £ ^ k ' (x " y) J h lV (h n (2-15) 

and 

Eui(k') = lov f sinfc' + cospf(1 — cosfe') (2-16) 

Notice we have dropped the important phase factor e lfco(5M from <?, for it plays no explicit role 
in the following analysis, since the limit N —> oo is taken before the limits L — >• oo and /3 — > 
oo. As consequence of fundamental properties of the Grassmann Gaussian integration, the 
decomposition of the covariance (|2.13p implies a decomposition of the field 

i>l, s = r£ ) + Y. E eiuPFe ^ e Sl (2-17) 



t,/3 



where fields with different scale labels or different label lo are independent, and the covariance 
of tp^ is while the covariance of tpu is <?w • Basically the label ui refers to either two 
branches of the dispersion relation. 

Let us now describe the perturbative expansion of the functional W( J, rj) defined in (|2.ip : for 
simplicity we shall consider only the case i] — 0. We can write: 

e w«o) = f F(# <0) f p^Wje-vW-^W+EJ^W = 

^ , ^ (2.18) 

= e -L(SE / p(^<0) e -V (D) W^)+8 <0) W i0 ,J) 



where, if we put x = (x x , . . . ,x 2n ), u = (oj 1 , . . . ,u) 2n ) and Vx,^ = 117=1 ^tuun IT»=«+i ^w.) 
the effective potential V^°\ip) can be represented as 

V(0) (^) = E E/ dx W S„(*)^ (2.19) 
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the kernels W^°2 n (x) being analytic functions of A and v near the origin; if \v\ < C\X\ and we 
put k = (ki, . . . , kan— i), their Fourier transforms satisfy, for any n > 1, the bounds, see §2.4 of 
[55], 

|wl°L(k)| < c ,n |A| max{1 '"- 1} (2.20) 

A similar representation can be written for the functional B^ a \ip-°, J), containing all terms 
which are at least of order one in the external fields, including those which are independent on 

The integration of the scales h < is done iteratively in the following way. Suppose that we 
have integrated the scale 0, — 1, — 2, obtaining 

e W(J,0) = e -L PEj J p Zi C ^< 3)e -V^\^^)+B^(^-^\J) (2 21) 

where, if we put C^k') -1 = J2h=h L $ AOO) Pzj.Cj is the Grassmann integration with propa- 
gator 

^^' ) (x-y) = ^^- £ ^" k(x ' y) -Xf\m ^ 
Z.j Zj /3L f— ' -iko + E u (k') 

V"'(^) is of the form 

vU) W = E E / ^<L(x)^ (2.23) 
n>l a 

and fiW (^— , J) contains all terms which are at least of order one in the external fields, including 
those which are independent on ip-i . For j = 0, Zq = 1 and the functional V*- -* and B^ are 
exactly those appearing in (|2.18[) . 

First of all, we define a localization operator in the following way: 



h d F x {$) +l 2 , j F 2 (,yZ~ j ij) + U ij F 4 (y/Z~iP) 



where 



(2.24) 



^ = E / ^V'x.^^oV'x^^ > ^ = ^E / dx ^>U,>* 



(2.25) 



x,w,srx,w,-srx,a;.-s 



and Thp XtU>a = Jrfke 4kx E u (k)ip£ >u>a (see definition ([2T6|)). Note that 

; 4 ,o = 2Aw(0) + O(A 2 ) 2 2 ,o = 2Av(0) + <3(A 2 ) Z 1)0 = 2Xv{2p F ) + 0(X 2 ) (2.26) 

and in writing (I2.24p the SU{2) spin symmetry has been used. In the case of local interactions, 
^(p) = 1. .Fi in (I2.25[) is called backward interaction while F 2 ,F4 are the forward interactions] 
the umklapp interaction, defined analogously as (|B.1[) below, is not present in £V^, as well as 
other terms quadratic in the fields. The reason is that the condition pp ^ 0, it says that such 
terms are vanishing for j smaller than a suitable constant (depending on pp and A), because they 
cannot satisfy the conservation of the momentum, so there is no need to localize them (more 
details are in ||24J. 
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Moreover, the local marginal operators associated with the densities (jTTTJ) are defined in the 
following way: 



r(«) 



£B (j) ( v r Z~ij,J) = J dx,l 
where 0( 1,Q ) are the small momentum transfer correlations, (see pag. 231 of |12| ) 



(2.27) 



(1,C) 

X 




V X 


LJ.S,S' 


n (i,sc) 

V X 


e.co.s 



(2.28) 



X,£*J, — S 



while 0( 2 <") are the large momentum transfer correlations, (see pag. 221 of (T2"] ) 

OJ,S 

(2,SC) = *V£ )Wl .V£,- W) -. 
s.uj.s 



(2.29) 



These definitions are such that the difference between — and —CV^' +£B^' is made 

of irrelevant terms. 

Note that the factor e " J£LJf>F in the definition of O x : comes from the fact that the two a £ 
operators in the definition (jTTTJ) of the triplet Cooper density are located in two different lattice 
sites (otherwise the density would vanish). Moreover, there is no local operator 0^ ,TC because 

S s s' u s^l^'V'x u> s ' = by anticommutation of the fermion fields. 

We then renormalize the integration measure, by moving to it part of the quadratic terms in 
the r.h.s. of JC3]), that is -z^L)' 1 s £ k [— jfc + ^(fc)]V>k iW , s Vw s ; equation (|2"72"Tj) takes 
the form: 

e w(,/,o) = e-Lpw+ti) J p _ ^ ^ W (<i)) e -v«(^^)+BO)(^^,^) (2 3Q) 

where is the remaining part of the effective interaction, i c (dtp-i) is the measure whose 
propagator is obtained by substituting in (|2.22[) Zj with 

Zj-xQl) = Zj[l + ZjCjik)- 1 } (2.31) 

and tj is a constant coming from the normalization of the measure. It is easy to sec that we can 
decompose the fermion field as ifi^-i) — ip^-i -1 ) + ipv\ so that 

P z,. u cM^ = Pz^c^Ad^-^Pz^jAd^) (2.32) 

where fj(k) (see eq. (2.90) of [35]) has the same support and scaling properties as /j(k). Hence, 
if we make the field rescaling ip — » [y r Zj~yj ^/Zj}^ and call V^(-y/ the new effective 
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potential, we can write the integral in the r.h.s. of (I2.30[) in the form 

J Pz^c^id^-^) J i 5 ZH! j :1 (# (,) )e 4U,( ^ ( - ,,)+fl ' ( ^ ,a ' J ' J) 

By performing the integration over ip^\ we finally get (|2 . 2 1 1) . with j — 1 in place of j. In order 
to analyze the result of this iterative procedure, we note that CV^(ip) can be written as 

d>0')(V>) - l j VjF v {i>) + 5 3 F a {^) + sijFiWO + 92 tj F 2 (i>) + <? 4 ,;W) (2.33) 

where vj = (y^/V^-iHn S j = {yf^il yj z i-\){ a i - z i) and 9i,j = (y^Zj/ 
i = 1,2,4, are called the running coupling constants (r.c.c.) on scale j. 

In Theorem (3.12) of [35] it is proved that the kernels of V^' and B"' are analytic as functions 
of the r.c.c, provided that they are small enough. One has then to analyze the flow of the r.c.c. 
(the beta function) as j — > —oo. We shall now summarize the results, following §4 and §5 of [24] 
with some improvement. 

2.3 The flow of the running coupling constants 

Define vector notations for the r.c.c, 

Vh = (Vl,h, «2,fc, «4,/u Vs t h> V Vl h) = {gi,h,92,h,94,,h,Sh,Vh) = (gh,Sh,Vh) ■ (2.34) 

The r.c.c. satisfy a set of recursive equations, which can be written in the form 

v QJ _i = A a v a j + /3^\vj;...,v ;\, v) (2.35) 

with A v — 7, A a = 1 for ci/i/. These equations have been already analyzed in [24], where it 
has been proved that, if A is real positive and small enough, then it is possible to choose v so 
that, fixed $ < 1, \v h \ < C\^ h , Vh < 0, and < g Ul < A(l + aAl^l)" 1 , for some a > 0, while 
the other r.c.c. stay bounded by CX and converge for h — > — oo. In this paper, in order to proof 
Borel summability of perturbation theory, we extend the proof to complex values of A, restricted 
to the set D e $ defined in (| 1 . 30p : this implies that we need an analysis a bit more precise of the 
flow equations (|2.35[) . 

To begin with, we put V\ = v and we suppose that the sequence {^}/t<i are known functions 
of A, analytic in D £t $, such that 

Kl < C\X\^ h , h<\ (2.36) 

and study the flow equations of the other variables. The idea is that this restricted flow has 
properties such that, by a fixed point argument, the sequence {i^h}h<i, satisfying the last equa- 
tion of (|2.35p . can be uniquely determined. Since this point can be treated in the same way as 
in spinless case (see §4.3 of [22]), we shall give for granted this result. Hence, from now on, we 
shall define Vj — (gi,j,g2,j,g4,j,8j) and we shall consider the restriction of (12.35)) to Vj. 

The next step is to extract from the functions fia the leading terms for j — > — oo. Observe 
that the propagator air of the single scale measure P 7 r-x, can be decomposed as 

^' ) (x) = ^^(x) + r«(x) (2.37) 

where <?q w is the Dirac propagator (with cutoff) and describes the leading asymptotic behavior 
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while the remainder rij satisfies, for any q > and < 1? < 1, the bound 

| r O)( x )|<2 7. is ■ (2-39) 

1 " v ;| ~ Z,- 1 + (7^|x|)« v ; 

Let us now call Z^j the values of Zj one would obtain by substituting V^ ^ with CV^ in (|2.18l) 

and by using for the single scale integrations the propagator (|2.3Tf) with r$ (x) = for any i > j. 
It can be proved by an inductive argument, see §4 of [35], that, if all the r.c.c. stay of order A, 



z j Z JM 



< Ce]^ (2.40) 



where 



max{|A|, max \g h \, max \5 h \} 



It is then convenient to decompose the functions pa as 

^\v 3 - ...,v ; A, v) = pWlyj, v ) + ^JHvf, ...,v ; A) (2.41) 

where ... , vo) is given by the sum of all trees containing only endpoints with r.c.c. Sh,gh, 

> h > j , modified so that the propagators gffi and the wave function renormalizations Zh, 
> h > j, are replaced by g^' u and Z^,h\ j3a contains the correction terms together with the 
remainder of the expansion. 

Lemma 2.2 

\^(v r ,...,v ;X)\<H^ r2 . ^ (2.42) 

As showed in [23], this lemma is basically a consequence of (|2.40[) and (|2.36l) . Therefore the 
leading term in (12.411) is /3a , that we further decompose as 

/3«)(v i ,...,v ) = ^(v i )+r a j(v i ,...,vo) (2.43) 

where /3a ^ (v) = (3a ( v , v )- We can write: 

^(v,-) = E & S( v i) + & ?> 2 ( v ^) ( 2 - 44 ) 

i=0,l 

where j(vj) is the contribution of order i in g\j, wile > 2 ( v j) is the contributions of all trees 
with at least two endpoints of type g\. The crucial property is the following lemma. 

Lemma 2.3 (partial vanishing of the beta function) 

|&S(v;)| < Cey* , i = 0,l (2.45) 

The above property was proven in §5.3 of [24], extending the proof for the spinless case in 
[37l [23] . and it will be reviewed in App. [C] Now, let us extract from Pa yVj) the second order 
contributions, which all belong to b^ >2 (yj); we get: 



/3«( Vj ) = -a^- + £ &S(vi) + ra,^) (2-46) 



i=0,l 
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with <zi = a > 0, a2 = a/2, 04 = as = 0, and, for some b\ > 0, 

1^(^)1 < hsjlgxj] 2 . (2.47) 
In the limit L, P = 00, if c^<7 = 5d,1> 

o r 1 f dk A>h) n \~(>h) n \ lo S7 /„ 

a - 2 ^ m oo h 7 (^f- 9d;+ (k) ®-- (k) = ^ ( 8) 

Let us now analyze in more detail the functions (v^-, Vq), which appear in (|2.43p . If we 
define, for j' > j + 1, 

(v,- , v ) = j8 W (v,- , .. , v s , v r , v„) - ( Vj - , ... , v j , vj , v ) (2 .49) 

we can decompose fa (vy, ...,Vo) in the following way: 



r ttJ (v„...,v ) = Yl ^(vj v ) (2.50) 

i'=3'+i 

Note that Da''' • Vo) is obtained from j3a\vj, — , Vo), by changing the values of the r.c.c. 
in the following way: the r.c.c. associated to endpoints of scales lower than j' are put equal to 
the corresponding r.c.c. of scale j; those of scale greater than j' are left unchanged; at least 
one of the r.c.c. v r ,j' is substituted with v r _y — u rj . By using the short memory property (see 
e.g. (4.31) of [35]), we can show that , if Ej is small enough, 

\D^(w j} v )| < b 3 e n -^'-^W r - v,| (2.51) 

for some 63 > 0. If we insert in the flow equation (|2.35[) the equations (|2.4ip . (|2.43[) , (|2.46|) . 
([2301 and use tlie bounds (j2~l2"1) . (POST) . (pl47|) and ([23T) , we get, if ej is small enough, 



~ vj\ < (a + 6 l£j )| 5 i J | 2 + (ceo + 62^)7* + 63^ £ l^'^Wf ~ (2.52) 

for some 02 > 0. The form of this bound implies that, in order to control the flow, it is sufficient 
to prove that gij goes to as j —00 so fast that |<7i j| is summable on j. Hence, we have to 
look more carefully to the flow equation of gij. By proceeding as before, we can write 

9i,3-i = 9i,j - a 9i,3 + h,j + n,j + fij (2.53) 



\ri,3\<b 2 sy\ Irkjl < b 1 s j \g 1 , j \ 2 , \njKbsSj £ y -*<f-i)\ Vj ,- v .\ (2.54) 

j'=j+i 

It is easy to show that, if £0 is small enough, there is a constant C4, such that, if gi^ € D SOi s 
and c 4 |j ||gi i0 | 2 < |^i,o| 2—T7 , V < 1, then, for j > j , 

\gi,a\/2<\g ld \<2\ gi , Q \, e 1 <2e (2.55) 

Hence we put jo = — (c4|gi.o| 1 ^ 2 ) 1 and suppose £0 so small that 

e j0 7* j0 <2c 5 |pi Jo |7 fj ' < l5i J0 | 3 (2-56) 
where we also used the fact that, since v(2pp) > 0, £q < C5 1 <7i,o| ? for some constant C5. 
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Lemma 2.4 If g± t o G D £0t $ and j > jo, t/ien, if Eq is small enough, 

|Vj_i - Vj-| < 2a| 5 ij| 2 + 2ce 7* j 



(2.57) 



Proof - We shall proceed by induction. By (|2.55[) . if Eo is small enough, ce + b 2 s 2 < (3/2)c£o 
and a + b\Ej < 3a/2; hence, (|2.57p is true for j = 0. Let us suppose that (|2.57l) is verified for 
j > h > 0. By (|235|) . it j >h> j Q , \gi,j\/\gi,h\ < 4; hence, by using (j2~52l and (j2~57j) . we get: 

|v h _i - v h | < (3/2)a\g Ul \ 2 + (3/2)ce ^ h + 
o 



max 



2a|3ijv| 2 + 2c£ 7- 



< |5i,fc| S 



(3/2)a + 64a6 3 eo^n7^ n 



+ 7* e<) 



(3/2)c + 4co 3 e ^n7~ f " 



n=0 



Hence, (|2.57p is verified also for j = h, if £q is small enough. 



The previous analysis implies that the flow is essentially trivial up to values of j of order 



l3i,o| 



-1/2 



(or even \gx,o\ v , < n < 1). If j < jo, we write (|2.53[) in the form 



and we define A* 



9i,j-i - 9i,j - a j 9i,j j 
and, for j < j 0l 
1 



a,4 = a — 



A 3 = 



Jo - .] 



j'=j+i 



9i. j = 



r i j + r i,j + r i j 



1 + A 3 9l,j Uo ~j) 



(2.58) 



(2.59) 



Lemma 2.5 There are constants c\,c 2 ,c^ such that, if gi.o € De ,s an d it £o * s small enough, 
then the following bounds are satisfied, for all j < jo . 



Ej < c 3 e 

|vj -Vj+i| < ci| ff i, J+ ip 

■ffij\<\ffi*\* /2 
- a| < c 2 \gi % 



\9i,j 
\a, 
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(2.60) 

(2.61) 
(2.62) 
(2.63) 



Proof - We shall proceed by induction. By using (|2.56l) . (12.571) and (|2.55l) . we see that the 
bounds (|2.60[) and (|2.6ip are satisfied for j = jo, if c 3 > 2, c\ > 3a and 4c£o < a - Moreover, 
9i,jo = Si.jo an d, by proceeding as in the proof of Lemma 12.41 and using (|2.56p . it is easy to prove 
that there is a constant c 2 , such that 



^Jo 



a\ < c 2 \gi 



Jo 



Hence, all the bounds are verified (for So small enough) for j = jo, if c\ > 3a, c 2 > c 2 and c 3 > 2. 
Suppose that they are verified for jo > j > h. 

The validity of (|2 .62[) for j = h — 1 follows from Prop. IA.21 which only rests on the bound 
(|2.63|) for j > h. On the other hand, f|2 . 62[) implies that, if £o is small enough, 2 _1 \gi j \ < \gij \ < 
2\gij\; hence, using (|2.59[) . we get, for j > h 



3i ,j 



9l,h 



< 4 



|1 + A h g ljQ (Jo - fa) | 
\l+Ajg hjo (j -j)\ 



(2.64) 
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Let us now define, as in App. [21 Aj = aj + i/3j, aj = $tAj, and suppose that 

2c2£o < a /2 



(2.65) 



so that, by (|2.55p . aj > a/2, \(3j\ < 2c2£o, \Aj\ < 3a/2, for j > h. By proceeding as in the 
proof of the bound (|A.8|) in App. [21 we get, if j > h and |Arg gi :0 \ < n — S, 6 > (so that 
|Arg 5lj0 | <n-S/2, see (233), 

I 1 + 9i,j (*j(jo ~ j)\ > jjj sin(<5/2)[l + |.ai j0 |a 3 - j)] 
and, if we put 1 + A,gi JO (jo — i) = 1 + a j9i,j Uo ~ i) + w ji we choose £q so that 



< 



6c 2 e \g ldo \(j -j) 



I 1 + 9i,j a jUo - j)\ sm(5/2)\g hl0 \{a/2){j a - j) 



12c 2 e I 



1(5/2) 



Then, by using (|2.64[) , we get 

9i, j 



9l,h 



< 24 l + (3a/2)|ff ljo |(jo-fe)| 
" sin(<5/2) l + (a/2)| 5l>JO |(io-j)| " * U ) 



(2.66) 



(2.67) 



for some constant C$, only depending on 6 and a. Moreover, since £h < C3£o, then c£o + bi£\ < 
2c£o and a + b\£j < 2a, if 

bic^e® < c , and bic^eo < a (2.68) 
Hence, by using the bounds (|2~52"|) . (j2~oT1) . (|2"3o) and ([2TF71) . we get 



|v/j_i -v h \ < 2a\gij l \ 2 + 2ce "f' 



-d(jo-h) 



\9i,jo\ 2 + cib 3 e h 1 
j=h+i 



^- h \j-h) max 

h<j'<] 



< \9i,hf 



2a + 2ce C 2 max 7 n€ n 2 + a£ h b 3 C 2 7 "'"n 

n>0 ^ — ' 



tfri 3 



n=0 



It follows that (|2.6ip is satisfied also for j = h, if 

00 

2a + 2ce Cj m&x-f-^n 2 + 2 Cl c 3 s b 3 C 2 V 7 -' 5n n 3 < c a 



(2.69) 



71=0 



Moreover, by using (|2.6ip and \gij\ < 2\gij\, we get, for some 64 > 0, only depending on a, 
under the condition (|2.65[) : 



£h-l < £0 + ^ l^" 1 _ V J'I - £ ° + fo 4Cl£ 

so that £^-i < c 3 eq, if 

1 + 64C1 < c 3 

The bound for ah-i — a can be done in the same way; it is easy to see that 



\a h -i-a\< 6ic 3 + 6 2 c3£ C|max7 "V + 2 Cl c 3 b 3 C 2 s V 7 

n > £ — * 

Hence, f|2.63p is verified for j — h — 1, if 

oo 

c 2 = 2cCj max7-"' 5 n 2 + 2ciC 3 & 3 Cf V 7~' 5n n 3 < c 2 

n>0 ^ — 4 



£o 



(2.70) 



(2.71) 



(2.72) 



n=0 
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The conditions ([2~BB1 . ([2351) . (12"1)9"1) . (|2T70)) and (|2~?2"j) can be all satisfied, by taking, for 

example, c\ — 4a, C3 = 1 + 4a&4 and C2 = max{c~2, C2}, if £0 is small enough. | 

Lemma 12.51 implies that, if gifi £ Z5 e < 5 , with e small enough (how small depending on 5), gij 
goes to 0, as j — > —00, and Ylj=h \di,j\ 2 — C<5 1 1 A | , uniformly in h. This is an easy consequence 
of the condition (|2.62p and the condition v(2pp) > 0; note that the power 3/2 in the r.h.s. of 
(I2.62p could be replaced 2 — 7], r\ > 0, but 2 is not allowed. The form of the flow (|2 . 35[) then 
implies also that g%j, 54, j and Sj converge, as j — > —00, to some limits 52,-00, 54,-00 and <5_oo 
of order A, such that 

52,-00 = 52,0 - i.9i,o + 0(|A| 3 / 2 ) = [20(0) - v(2p F )} A + 0(|A| 3 / 2 ) (2.73) 

54,-oo = 54,o + 0(A 2 ) = 2Xv(0) + 0{\ 2 ) 
S-00 = O(A) (2.74) 

Let us now suppose that A is a (small) positive number; the previous bounds imply that g± t j > 0, 
for any j < 0. The following Lemma will allow to control the logarithmic corrections to the power 
law fall-off of the correlations. 



Lemma 2.6 There are four sequences Wi^h, <5i.h, i = 1, 2, h < jo, such that 



30 ^ 

^2gi,j = (1 + wi,h)- log[l + ag ljjo (j - h)] + S ljh 



3 = h 



[92,j ~ 52,-oo] 

3=h 



(1 + w 2 ,h)— log[l + agi ijo (j - h)] + 5 2 ,h 



with 



M<cx, im^a 1 / 2 

CX 



[1 + ag ltjo (J - h)\ log[l + ag ly3o (j - h)} 



Proof - Let us put 50 = 5i,j j and a(s) the function of s > 0, such that a(s) 
n < s < n + 1 . Then, by using (|239| . (l2~62)) and ([2~63)) . it is easy to see that 



(2.75) 

(2.76) 

(2.77) 
(2.78) 



a jo—ni if 



3o 



j=h 



< CX 1 ' 2 , I n = ds- 

JQ 1 



5o 



50 Jo dt a(t) 

On the other hand, (|2.63l) also implies that a(s) — a + Xr(s), with \r(s)\ < C; hence 

5() 



(2.79) 



In = 



ds- 



X 



9o Jpdtrjt) 



1 + 5oas J [1 + g Q J Q dt a(t)} [1 + g as] 



implying that 



/„ log(l + ag n) 

a 



< 



4CA f a9on 



dx 



4CA , „ 
< ~r log(l + ag n) 



Hence there is a constant w n such that /„ = (1/a + w n ) log(l + agon), with \w n \ < CX; this 
bound, together with the bound in (I2.79[) , proves (|2.77l) for i = 1. To prove (|2.78p . note that 



\In+l — In\ < 



n + 1 



ds 



go = £ 

l + 5ofs a 



log 1 



Q5o 
2 + agon 
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'n+1 — J-n 



(l/a + w„+i)log 1 + 



ago 



1 + ag n 



+ (w n+ i - w n ) log(l + ag n) 



so that, if A is small enough, 



|io n+ l - w n \ log(l + ag n) < + CXj log M + — a, " ! 



< 



45o 



+ agon J 1 + ag^n 
To prove (|2~77ji and ([2TT8")) for i = 2, note that, by (pOSI) and Lemma [231 if j < j , 



ff2j - 52,-oo 



9i,h 



h=—oo 



\ + O(A) 



b -| (1 + ag s) ; 



2+ 0(A) 



Hence, the proof of (|2.77l) is almost equal to the previous one, while the proof of (|2.78[) needs a 
slightly different algebra; we omit the details. | 



2.4 The flow of renormalization constants 

The renormalization constant of the free measure satisfies 

%i = 1 + ^ %,Si, g , So) + W ] (%■; .., v ;X) ; 



(2.80) 



while the renormalization constants of the densities, for a = C, Si, SC, TCi and i = 1,2, satisfy 
the equations 



Z 



(i,a) 
2-1 



^) = 1 + ® ^2 - - , 9o, So) + ^ (vj ; .., v ; A) 



(2) 



(2.81) 



In these two formulas, by definition, the /3{ functions, with £ = z or (i, a), are given by a sum 
of multiscale graphs, containing only vertices with r.c.c. gh,Sh, > ft. > j, modified so that 
the propagators and the renormalization constants Zh, , > h > j, are replaced by 

5d*L' ■^i £> ' ) ' Z^' 1 ^ (the definition of z\^' l,a " 1 is analogue to the one of Z^); the /3p^ functions 
contain the correction terms together the remainder of the expansion. Note that, by definition, 



the constants z\ D> are exactly those generated by (j2~80| and (j2~4T1) with /3^ j = f3^> = 0. Note 



(2) _ 5(2) 



also that |/3p } | < C^V^', while |/3« j q) | < C%7 
By using (I2.80[) and (|2.81|) . we can write 



Z 2-i 
^2-1 



+ ^l Q )(52,^) + /3^ 1 , Q) (^;..,«o;A) 



(2) 



(2.82) 



with |/3 z ,(i, Q) | < Cv ^K If we define P ( ^ {l a) {g,5) the value of ^( 1>a) (Sj, 5*0, 5 ) at 

(g*i,8i) = (g,5), j < i < and /S^'q „\ (if, 5) the sum of its terms of order and 1 in gi t h, it 
turns out that 



\& j) (^,Sj)\ < C[max{| 5l || 52 |, | 54 |, |*|}]V" , if « = C 



(2.83) 
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This bound, as crucial as the analogous bound (|2.45p . has been proved in [33]; the proof will 
be sketched in App. [Cj The bound (|2.83l) . together with Ylk=j I *7i , | 2 < C|A| and the fact that 



Z 



(l,Si) _ 7 (1,C) 



by the SU(2) spin symmetry, imply that 



<C\e 2 A, a = C,S z 



(2.84) 



Regarding the flow of the other renormalization constants, we can write 



(*) 



(2.85) 



f z) 

where Z ■ = Zj and, by definition, 



r\t = hm 77^- := log 1 + /3 4 (0j) (52,-00, 54,-oo, <5-oo; -5 52 -00, 54 ,-00, £-00) 

j->-oo L 



(2.86) 



Note that the exponents rj t are functions of U-00 only, an observation which will play a crucial 
role in the following. Moreover, by an explicit first order calculation, we see that 



' {2m F )- 1 g 2 ,- O0 + 0{\ 2 ) t = (2, O, (2, Si) 
//, = ^ -(2 7 r t ; i ,)- 1 52 ,_ 00 + 0(A 2 ) t = (2, SC), (2, TQ) 
0(A 2 ) otherwise 



(2.87) 



while 



Z 



h-l 



7(2,0) 
7 (2,C) 
7 (2,5.) 
7 (2,S<) 
A(2,SC) 
£(2,SC7) 
7 (2,TC 4 ) 
7 (2,TC ( ) 



1 + 0(ft, fc A) + r^> , i = z, (1, a), a ^ TC 4 



(2,C) 



1 - a.9i,/i + ^ (52,A - .92,-00) + 0(5i, h A) + r}; 



1 + |(fl3,h - .92,-00) + OQfi jh A) + rf Sl) 



1 - |si,/t - ^(92,h - 52,-00) + 0(5i,fcA) + ^ 2,SC) 



1 + ^9i,h - ^(g2,h - .92,-00) + 0(g hh X) + rf' TCl) 



where a and 51^ are defined as in (I2.48P and (|2.59p . respectively, and X)°=~oo l r l^l — C j A| s 
Let us define: 



CO 
9t = 



log(l + ag lfi \h\) 
Hence, by using (|2~75|) . ([2~7S]) and (|2~771) . we get 



(2.89) 



19/ 



(ft) 



|g£ fc) |<CA, t = z,(l,a),a^TCi 
ka|<CA, t=(2,a) 



(2.90) 



where the constants £ Q are those of (|1.20[) . 
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2.5 Proof of Theorem fTTTI 

In order to prove the representation (|1.19l) of the density correlations J7 Q (x), we can proceed 
exactly as in §5 of [35], where a similar problem was treated in all details; hence we shall only 
describe the result. We can write for £7 Q (x — y) a convergent tree expansion, whose trees have 
two endpoints associated to density operators (special endpoints) and an arbitrary number of 
interaction endpoints (normal endpoints). As one could expect, f2 Q (x — y) behaves, as |x — y| — > 
oo, as the function f2 a (x — y) calculated by taking only the trees with two special endpoints of 
scale h < and no normal endpoints. This function can be obtained by the following procedure. 
Let us consider the expression 

(0^;0^)l + {0^;0^)l 

where the operators Ox'"', i = 1,2, are those defined in (I2.28P and (12.291) and ( • )d i s the 
truncated expectation evaluated with covariance Y2h Z^ 1 g^ LJ . Hence, this expression can be 
written as a sum of terms, each one proportional to Z7 1 (^^(x — y) Z^g^ ^,(x — y), for some 
values of h, h', lu and u/. tt a (x. — y) is obtained by multiplying each one of these terms by 
i z hvt'] 2 (hVh' = max{h,ti}), if it appears in the calculation of (Ox Oy' Q )p, otherwise by 
[^ftvft'] 2 - Let us consider first the case a = C; we have 

Q c (x) = (1 - c) (x) + cos(2p F x)n (2 ' C) (x.) 

f^)(x) = 2]T]T [ -§f^9^ M<?£2(x) (2.91) 

w h,h' h h ' 

n^( X ) = 4 £ [ ^f^9<£l (*)<#>(*) 

h,h' 

Let us now observe that, for any N > 0, |sj^(x)| < C N j h [l + (7 /l |x|) JV ]- 1 . Hence, if |x| > 1, 
in the previous sums the main contribution is given by the terms with \h\ and \h'\ of the same 
size as log 7 |x|, so that one expects that the asymptotic behavior of Q( l ' C )(x), i — 1,2, is the 
same of the the function f2^' C ^(x), obtained by the substitutions of r )~ h and with |x| in 
the asymptotic expressions of the renormalization constants, given by (|2.85p and (|2.88[) . that is 

Ui^L^|x| 2 ^-^ l + /(A)log|x| V C - ' (2.92) 

where the coefficients q\ are defined as in f|2 . 89[) . h x — ini{h : 7 ,1 |x| > 1}, and, by (|2.26|) . 
(j2~59l and Lemma 

a^^2_Xv(2p I ) + , 
log 7 itvf 

In order to justify the substitution (I2.92[) . let us put r\i = 2(iji^c — f] z ) and <7i(x) any continuous 
interpolation between 2[qf^ — qi h ^] and 2[q^ ^ — qi hx Note that, thanks to the bounds 
(|2.77[) and (|2.78p . <Zi(x) is a bounded function of order A, defined up to fluctuations bounded, 
for |x| > 1, by CA[L(x) logL(x)] -1 , with L(x) = 1 + /(A) log |x|; hence, its precise definition 
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modifies the following expressions only for a factor 1 + O(A). Let us now note that 



|fiC.O(x) - ^< c )(x)| < CW|x|"*- 2 [l + /(A) log |x|]*« £ - 



7 x 



7 fc '|x| 



(7^x1)^(7^1x1)^ 



(yiV/i'| x |)277i,c 



£(|x|) 



L( 7 W) 



/i./i 

,_ ,/,( fc ') - ..... o_(kVh') 



+ (7*|x|)"l + ( 7 fc'|x|)" 



£(|x|) 



L( 7 lft'l) 



^(|x|) 
L(7l hvh 'l) 



-29];: 



(2.94) 



where 



L(t) = l + /(A)logt, c fe = L( 7 l^ h, /^\ c h = L( 7 |,(| )-^/^ i ' C) 

By ([27gg]l . (1277511 and (127751) . c ft = f + 0(A 1 / 2 ) and 4 = f + 0(A 1 / 2 ). On the other hand, if r > 
and t 7^ 0, 

|r* - 1| < |ilogr|(r* + r~ 4 ) 



and, if g 7^ 0, 



i(|x|) 



L(j\ h \) 

These two bounds, together with the bound 



|/(A)log( 7 k |x|)| + |/(A)log( 7 *|x|)|l«l+' 



( 7 'V) Q ilog( 7 'V)| /3 



+ ( 7 ' ! T) Ar 



5: CN,a,q 



h — — oc 

valid for any j3, r > 0, a > and N > a, imply that 

|f!^ c >(x) - ^ c )(x)| < CWA^Ixl*-^! + /(A) log |x|]*<* 



(2.95) 



By the remark after (|2.59[) . the factor A 1 / 2 can be improved up to A 1- " 5 , i? < 1. 

By proceeding as in §5 of [35], it is possible to prove that a bound of this type is satisfied also 
from the sum over all the other trees. Hence, in order to complete the proof of (|1.19|) in the case 
a = C, we have only to calculate ^^'^(x) and Q,( 2 ' C '(x). By using l|2.84[) . we see that rj\ t c = 7 h 
and q[ h l, = q^, so that, if we define X c = 1 - m,c - Vz an d Cc( x ) = 2[g 2 ,c(x) - &(x)], we get 



n( 1 ' C )(x) = 2^ ffD ^(x) 5D , w (x) 
17( 2 ' C )(x) = 4|x| 2 ( 1 -^)[l + /(A)log|x|]^W 5D , + (x).ga_(x) 
where <?d,w(x) = J2h=-oo 9d u( x )- ^ n tnc °t ner hand, it is easy to see that, for any N > 2 



(2.96) 



1 



1 



+ o(\*\- N ) 



2ir vf%o + iu>x 

It follows that, up to terms that we put in the "remainder" Rc(x), 



^> c >(x) = 4^o(x) , ^>(x) = ^^ 



(2.97) 



where the functions f2o( x ) and L(x) arc defined as in Thcorcm ll.il Hence, by using (|2.87p and 
(|2.93[) . we get f) 1 .19(1 for a = C, together with the fact that Cc( x ) = —3/2 + O(A), in agreement 
with (fl~20| . and 2X C = 2- b\ + 0(A 2 ), in agreement with (fTT22|) . Note also that, in Theorem 
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11.11 we have modified the function /(A) by erasing the terms of order greater than 1 in A; the 
only effect of this modification is a change of the function Rc^), which does not change its 
bound. 

The proof of (|1.19[) in the other cases is done in the same way. In particular, in the case a — Si 
we have to use again the bound (|2.84l) , while the fact that there is no oscillating contribution to 
the leading term of firCi is due to the fact there is no local marginal term which can produce 
it, by the remark after 92.291 

Finally, the proof of the scaling relations (|1.22[) follows from the important fact, proved in 
93.61 that they are the same as those of the effective model. Hence they follow from the explicit 
calculations of 93.51 



First consider the free energy (|1.4p . We can decompose it as E(X) = Y^h=-oo Eh(ty, where 
Eh(\) is the contribution of the trees whose root has scale h, hence depends only on the running 
couplings Vj with scale j > h. The tree expansion implies that there exists £q, such that, if 



then \Eh\ < C2"f 2h £o, with C2 independent of h. The analysis of 92.31 implies that, given 5 G 
(0, 7r/2), there exists e such that, if A e D e> $, the condition (12.98[) is verified uniformly in h; 
then it is easy to see that E(X) is analytic in D e ,s and continuous in its closure. The domain 
of analyticity of Eh(X) is in fact larger; the form of the beta function immediately implies that 
there exist two constants C3 and c such that Ao < C3 1 A| and, if A/j < £0, then Xh-i < Aft + cA^; 
hence, if c 3 |A| < min{e /2, l/[4c(|/i| + 1)]}, then, if j > h and Xj < 2X , 



It follows that Eh{X) is analytic in the set (11.311) . with cq — C3 1 min{eo/2, l/(4c)}, and that 
\Eh(X)\ < ci7 2/l , with ci = c 2 £o; hence E(X) satisfies (ll.32j) with n = 21og7. 

Let us now consider the 2-point Schwinger function 52(x). By using the tree expansion 
(similar to that written in [20] for the infrared part of the spinless continuous Fermi gas), we 
can write ^(x) = J2h=-oo ^^(x), where S2,/t(x) is the contribution of the trees whose root has 
scale h. By proceeding as in §6 of [20], we can prove that, if (|2.98l) is verified (possibly with a 
smaller eq), then, for any N > 0, 



with c/v independent of h. Hence, if we define h x so that 7' ix |x| E [1,7), then, if /i x > h 



2.6 Proof of Theorem [Q] 



X h = max \vj\ < £ 



(2.98) 



X3-1 ^~ X o + b'lcA - < A (l + 4|i|cA ) < 2Ao 




(2.99) 




(2.100) 



h=h+l 



h—h x 



and a similar bound holds for h x < h so that 




and we can proceed as in free energy case, so proving (|1.32[) for S^x), with c\ = c s (l + |x|) -1 / 4 
and = I (this value could be improved up to any value smaller than 1, at the price of 
lowering £q down to 0). 
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The previous argument can be extended to the generic 2n-point Schwinger function, by using 
the same strategy used in §2.3 of |38j to analyze the corresponding tree expansion in the case 
of the Thirring model. The proof of the Theorem in the case of the density correlations is very 
similar to that of the 2-point function case, if one uses the description of the tree expansion given 
in §5 of [35]. We shall not give any further detail; the idea at the base of the proof is always 
that, in the tree expansion of the correlations at fixed space-time points, the contribution of the 
trees with the root at scale h must decrease exponentially with the distance from some fixed 
scale depending on the space-time points. 



3 RG Analysis of the Effective Model 
3.1 Introduction 

We introduce an effective model (related to the g-ology model introduced in |12) but differing 
because the interaction is time-dependent), describing fermions with linear dispersion relation 
and a non-local interaction; such model will be used to prove the crucial bounds (I2.45[) and (|2.83l) 
(on which the previous analysis is based) and for the proof of the Luttinger liquid relations (|1.22[) . 
together with Theorems ll.2l and ll.3l The model is expressed in terms of the following Grassmann 
integral: 




where x g A and A is a square subset of K 2 of size 7 , say 7 /2 < |A| < 7 , Pz(di/j l ' N ) is the 
fcrmionic measure with propagator 

k 

where Z > and S are two parameters, to be fixed later, vf is defined as in Theorem 11.11 and 
Xi,N(t) is a cut-off function depending on a small positive parameter e, nonvanishing for all k 
and reducing, as e — > 0, to a compact support function equal to 1 for 7' < |k| < j N+1 and 
vanishing for |k| < r ) 1 ^ 1 or |k| > •y N+1 (its precise definition can be found in (21) [22 ]); 7' is 
the infrared cut-off and 7^ is the ultraviolet cut-off. The limit N — > 00, followed from the limit 
/ — > —00, will be called the limit of removed cut-offs. The interaction is 

V(il>)=g 1 , x V 1 ,xW+g\\V\\{i>)+gxV±(fl)) + g 4 y^) (3.3) 

with 



vv(V0 = 


UJ,S 


f <ixdy/i(x 


- y)^x, w 


lS V>x, W 






W) = 


\y J 


f dxdy/i(x 


- yfe 








VxW = 


2^ J 


f dxdyh(x 


- y)< u 




%b + 




W) = 




f dndyh(x 


- y)Vi w 




,sYy,u>,- 





(3.4) 
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where h(x — y) is a rotational invariant potential, of the form 

Mx-y) = -^EMP) e4p(x ~ y) . ( 3 - 5 ) 
p 

with |/i(p)| < Ce~^ p \, for some constants C, fi, and h(0) = 1. 

The model with g±,± = is invariant under the global phase transformation ift^ s — > 
e ±ia^, s% p± ^ ^ with the constant phase ct W)S depending both on u> and s. However, if gt t ±_ 7^ 0, 
the model is only invariant under the transformation ip^ w s — > e ±la "Vx o> s with the phase inde- 
pendent of s. 

The removal of the ultraviolet cut-off is controlled by an easy extension of the analysis given 
in §2 of [36] or in §3 of [39] for a spinless model with interaction AVn (ip); the presence of the other 
terms (including the (73-interaction considered in App. [B]) produces more lengthy expressions 
but introduces no extra difficulties. The crucial idea is to use an improvement respect to the 
power counting bounds due to the non- locality of the interaction, and to use that the " fermionic 
boubble" (see (2.39) of [5B] or (3.17) of |28| ) is exactly vanishing. 

Regarding the removal of the infrared cut-off, we perform a multiscale analysis very similar 
to the one given in <J21 that we shall sketch in App. [B] and App. [D] (we will refer to §4 of [2~4] 
for more details). It turns out that the infrared cut-off cannot be removed by this technique for 
all values of the couplings, but we are able to consider only two situations leading to a bounded 
flow: 

1. the case gi ^ = 

2. the case g\\ = g± — gi.± and gi t ± > 

In the first case, when the ultraviolet and infrared cut-offs are removed, the model becomes 
exactly solvable, a property related to the invariance under the local phase transformation 
w s — * e±iQic,I * J ' s ' ( /'x,w,s- Indeed, as we will see in ^3.21 - ^3.51 the functional integrals generating 
the correlations can be exactly computed, up to corrections which are proved to be vanishing in 
the removed cutoffs limit. This will allow us to prove (I2.45[) and that the exponents of this model 
and the Hubbard model (also analyzed via functional integrals) are the same, if the parameters 
of the effective model are suitable chosen (see ^3.61 below), so that the universal relations (jl.22l) 
follow. In the second case the model is not solvable, but still some correlations can be exactly 
computed, see Sj5] and this, again via a fine tuning of the effective model parameters, allows us 
to prove the relation Q1.24[l . 



3.2 Ward Identities in the g lt ± = case 

In the gi t ± — case we can derive a set of Ward Identities (WI). If we make in the generating 
functional (|3.1|) the gauge transformation w s — > e'° w ifj US , we obtain, in the limit of 
removed cutoffs, the functional Ward identity (WI): 



d m(p )^^-i 3 _,Cp)E^(p) 



dW (J, 77) 
3? 

L/t -'p, — a.r 



AtiZ 2 c 



(3.6) 



where -D M (p) = — ipo + coop, 



1 f dk 



dW 



dW 



drj 



k+p,/x,s 
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and 

^"(p) = [<W (8 a ,-ig± + S 8 ,ig\\) + S u -iS s 

Summing (|3.6[) over s we obtain the charge Ward identity: 



-104 



Hp) 

Aire 



(3.7) 



£> M (p)-i/4(p)r>- M (p) 



E 



2i/„(p)£>_^(p) 



<9W 



E 



-£>-m(p) 



AttZ 2 c 



+ S PjMjS (J,r?) 



(3.8) 



with 



i/ 4 (p) = 34/i(p)/(4ttc), 2zy p (p) = (.g,, + g±)h(p)/(4nrc) 
Multiplying (|3.6p by s and summing over s we obtain the spin Ward identity: 



£>„(p) + !/ 4 (p)-D-^(p) 



E s — 



2z/ (7 (p)£>_ M (p)^s 



P,M>s 



<9J, 



P,-M,s 



AttZ 2 , 



+ B P)M>a (J,rj) 



(3.9) 



with 

The proof of {3H and 



2^(p) = (ff|| - flx)/i(p)/(47rc) (3.10) 

with J = is essentially identical the one in §2 of [35J or in §3 
of [35]) while the presence of the term linear in J in the r.h.s. is explained in Sec. IV B of [ID] 
(see also App. A of [23); hence, we will omit here the details. 



Note the presence, in the WI p.8[) and (|3.9p . of the 



'ai "pi 



z/4 terms, which are called 



anomalies] they appear as a consequence of the breaking of local symmetries in the functional 
integral (|3.ip . due the momentum cut-off xi,nQ*) in the fermionic integration. This symmetry 
breaking produces extra terms in the WI which do not vanish when the cut-offs are removed. 
Note also that such anomalies are linear in the coupling. Such a property is called anomaly non 
renormalization, and is crucially related to the non locality of the interaction [36 ; in presence of 
local interactions, like in the massless Thirring model, it can be violated |38] . Another important 
point to be stressed is that (|3.8p is true also when g± ± > 0, while (|3.9[) is not; this remark will 
be used in the proof of Theorem 11.21 

By some other simple algebra we obtain from p.l3[) . (|3.14[) the identity 



dW 



MZ,Jp) + s> S M° 
^ 2 

f.l.S 



.(P) 



where, if 7 — p,a, and setting v± iP 
~D-„(v) 



J- 



4ttZ 2 < 



(3.11) 



-2/4,, 



1/4, 





D- P (p) - 1/4,7 (p)-D^(p) 


<W + 


2i/ 7 (p)£> p (p) 


V.-p 


D+(p) - J/4,-y£>- (P) 

w 7 ,+(p)<W + w i-+(p 


D_(p) - ut t7 (p)D + (p) 

, U 7,-(p)<V,M + 


-4^(p) jD+ (p) jD _(p) 

w-t ,-(p)V -11 



(3.12) 



"*«7,+ (p)(P0 + «>U 7 (p)cpi) 



for 



%,p(p) = 5 



1-^4, 7 (P) , 1 



V. 



■yy 

,+(p) 
1 



/'- 



V4 , 



-(P) 



/'-- 



«7,+(P) U 7,-(P) 



(3.13) 
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r,/*(p) = ( : ~ ^4, 7 (p)) - 4^( p ) , U7 (p) = u 7) _(p)/« 7i+ (p) 



(3.14) 



By doing suitable functional derivatives of the functional WI 



we can get many WI 



between the correlation functions. For example, if we make two derivatives w.r.t w s and 
a; s m both sides of (|3.8[) . we sum over /i and we put 77 = J = 0, we get the following charge 
vertex Ward identity: 



-Wo 



where 



1 - 2^ p (p) - 1/4 (p) 



(p;p 



1 - 2^ p (p) + 1/4 (p) (/,: ,;p:p k) 



G 



= ^[G 2 ;^(k)~G 2;w , s ( P + k)] 



<9 3 W 



s (p;p + k) = 
Giio,,.(p;p + k) = 2^A*r^ r^f — 



G2.U 



d 2 W 



(3.15) 

(3.16) 
(3.17) 
(3.18) 



In a similar way, the functional WI (|3.11l) can be used to obtain a closed expression for the 
correlations of the density operator p x ,w,s = V'x w sV'x w s- I n fact, if we take p. Ill) with 7/ = 
and we perform a derivative w.r.t. J_ 



we obtain: 



(Pp,. 

which implies that 



'P-p,< 



-X>- U (p) 



Mp) + *'^C.Jp) 

4vrZ 2 c 2 



where 



(/>x,w' 



1 r 



/T " 2 
1 

4ttZ 2 c 



(3.19) 



(3.20) 



^P ,.p(x-y) P0+cV 



(27T) 



A,(p) 



3.3 Schwinger-Dyson and Closed equations. 



By substituting the Ward Identities found in the previous section in the Schwinger-Dyson equa- 
tions, one get a set of closed equations, up to corrections which are vanishing in the limit of 
removed cut-offs; this is due to the non locality of the interaction and the proof is essentially 
identical to the one in [3(5] or in §4 of [3J5] for the spinless case. The presence of the spin makes 
however the resulting closed equations much more complex, and new properties emerge from 
their solution, like the spin-charge separation phenomenon. 

Given any F(ijj) which is a power series in the field, the Wick Theorem says that, if < ■ > 
denotes the expectation w.r.t. the free measure, 



(^,^)>o=5D^(k)(|^) . 



It follows that 



de 



vv 



k,u,s 



ffD,"( k )( 



(3.21) 



(3.22) 
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Hence, by using (|3.3D and (13.71) . we find, in the removed cutoffs limit: 

de w 



^(k)^(0,,) = ^ k _e w ^) 



dp 

7^ W 



(p) 



2 C W 



9 2 e 



(3.23) 



-(0,r?) 



that we shall call the Schwinger- Dyson equation. Let us now take the WI in the form (|3.11j) . 
with J = and \j! = —u>', derive it w.r.t. ?7p~ +k u s and insert the resulting expression in (|3.23[) : 
we obtain the following closed equation 



A,(k) 



de 



1 



Z 

UJ.S 



w 



E 

fx . t 



d 2 e w 



(2tt) 4 ~^ s * vp; 



lq,fj,,t" '/k+p,a;,s 'k+p.w.s 'q+p,Ai,i 

where, if is defined as in (|3.12[) , 

3.4 The two point function 

By using (|3.24[) . we easily get: 

(t^ijy,^) = <W<S»,*"SUx - y) 
where <S w (x) is the solution of the equation: 

(dM (x) - Fr w ,+ (x)5 w (x) = ~<J(x) , 
with d w = 9 Xo + iojcd Xl . The solution of (|3.25p is: 

5 u (x) = e A ^ x l°' 5w (x) , ffu (x) = -±- ^- , 

2nZ cxo + iwx 

having defined such that <9*A^(x|z) = Ft u s ( x ): 



Af(x|z) 



(2tt) 2 £> w (k) 



-F* (-k) = A*(x|z)+ S A*(x| Z ) 



for 



A»=^(-p) 



(3.24) 



(3.25) 



(3.26) 



(3.27) 



(-P) , g n , , K-o, e (-p) 
A,(p) + 2 ftl pj D u (p) 

Act(p) = ^ p) ^(p1 T* ( - p) D U P) 

In order to evaluate the asymptotic behavior of Af (x|0), we need to study functions of the type 



d 2 p 



o(p) 



e -»p x _ 1 



(27r) 2 (po + «wcpi)[u+(p)p - isujv-(p)cpi 



(3.28) 
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where o(p) and v s (p) > are even smooth functions of fast decrease. It is easy to show that 

o(0) 



Wx) = -^/^(x) +A + 0(l/|x|) 



«+(0) 

where A is a real constant and, if v = u_(0)/tL|_(0), 



(3.29) 



+1 dpi f + °° dp Q e-iboxo+pisi/c) _ i 
l-i ( 27rc ) J-oo ( 27r ) (Po + iupi)(po - ieujvpi) 
One can see that, if v > 0, v ^ 1 and x^O, 

1 



2ttc(1 + ev) 



[F(x ,ujxi/c) + eF(vx 0l —euxt/c)] 



where 

F(x ,xi] 



1 *i L-Pi(ko 



t+isgn(a;o)a:i) 



Pi 



In \z\ + iArg(sgn(x )2:) + B + 0{l/z) 



where z = xq + zxi, -B is a real constant and |Arg(z)| < it. Since 

Arg(sgn(a;o)2;) = Arg(z) - (x )sgn(xi) tt 

the function F(x) (considered only for |x| > 1) is discontinuous at xq = 0, while 7 WiE (x) is 
continuous. We can then write 



1 



2ttc(1 + ev) 



[log(xo + iujxi/c) + elog(wa;o — iewxi/c)] + C + 0(l/|x|) (3.30) 



where C is again a real constant. By using (|3. 12|) . Q3.28p . (|3.29[) and (I3.30p . one can easily check 
that 



A 7 (x|0) 



In (v y x a + (xi/c) ) ! — ; — In 



for 



Aire "' y 1 

+ q + o(i/|x|) 

7,6 v ~t,+ + £v -r,- 
2<7 7 w 7>e + <74 j7 w 7:£ 47re 



47TC 



w 7 x + iujxi/c (3.31) 



£5 7 



i> 7i +u 7j _ 

K - - CT 7,+) S 



2^ 7j+ j/ 7 ,_ 



1 - 



v j,+ ~ £v -y,- 

where are real constants and u 7 = w 7i +(0)/w 7) _(0) (and 34, p — gn while g4 jCT = —54). 
By using (|3~T5)) and ([3~Ti| . 



7: + 



C7 + ^ 



7,+ 



47TC V 7i + 1> 7: _ 



47TC 



= 



Aire 



1 ~ |Kv+ +^7 ,-) 2 _ ^7 g 7 - + g 7,+ 

2 4ttc 



(3.32) 



(3.33) 



2w 7]+ u 7j _ 

Note that this expression is continuous in u 7 = 1, as one expects, and that, at least at small 
coupling, ?7 7 > 0. 
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By using (|3.26[) and (13.271) . we finally get 



a M _ 1 (c 2 ^xg + afl-Wa^g + x 2 )^/ 2 

where C is a real constant O(g) and z 1 / 2 = \z\ 1 / 2 e lAr9 ^/ 2 . Note that the leading term is well 
defined and continuous at any x^O. 

Note also that, if (74 = 0, v p = v a = 1 and r\ p = r\ a = r)/2, so that 

= 1 (c 2 , 2 +, 2 )-" /2 e c +0 (i/|x|) (3 35) 

27TZ CXo + ZWIl 

If we also put gj^ = and £/|| = A, we get for r\ the value found for the regularized Thirring 
model, that is 77 = 2t 2 /(1 - t 2 ), with r = A/(4ttc); see eq. (4.21) of [39]. 



3.5 The four point functions and the densities correlations 

We want to calculate the truncated correlations (O^Oy' > ) T of the local quadratic operators O x , 

t = (l,a) or (2, a), defined as the analogous operators of the Hubbard model, see f|2.28[) and 

(I2.29|) ; note that pf has no special meaning in the effective model, but it is left there since we 

want to compare the correlations in the two models. 

Our UV regularization implies that (O x ) = for any t; hence we can make the calculation 

very simply, by using the explicit expressions of the four points functions which follow from the 

closed equation (|3.24[) and then evaluating them so that the two coordinates corresponding to 

each O^' operator coincide, if this is meaningful. This works for all values oft, except (1, C) and 

(1, S3), where there is a singularity, related to the fact that the operators p Xl u, s = sV'x u> s 

are not well defined in the limit N — ¥ 00, because of the singularity of the free propagator at 

x = 0. However, in these cases we can use directly the WI (|3 . 19[) for the density correlations, 

( 1 CI ( 1 s 1 

which allows us to calculate correctly, in the limit N — > 00, the correlations of O x ' and O x ' , 
by using ([3T2D]) , (|3~i"2")l and the equations (|3T55|) . ([5^9]) . (|53D]l . We get, for |x| > 1, 



1 - v. 



8n 2 c 2 Z 2 
2 



' ~ %k 2 c 2 Z 2 



y 7 ,+ — t> 7 ,_ (v-yXo + iojxi/c) 2 t> 7 ,+ + w 7i _ {v^xq — iujxi/c) 2 
w 7 ,+ 1 w i,- 1 



v 7 ,+ — u T) _ (w 7 a;o + iuJXi/c) 2 v 7j + + t! 7 ,_ {v-yXQ — iujxi/c) 2 



the corrections being of order l/|x| 3 . This implies that, for |x| > 1, 



while (O^ ' a; Ox T is obtained from this expression, by replacing 1/4 with —1/4 and v p with 
i/cr (hence also v p , v Pi + and u P! _ with v a , u CT) + and i>o-,-)- 

One can see, see Lemma 4.1 of [H], that the same result could be obtained starting from the 
four point function, if we take the limit e — > of the expression obtained by the substitution of 
the density operator p x .w.s with the regularization 



Px,^ s = J duS E (u - x)VXs<X 



where 5 e (x) is a smooth approximation of the delta function, rotational invariant (in agreement 
with our UV regularization), whose support does not contain the point x = 0. 



33 



In order to calculate the other correlations, we first note that the only four points functions 
different from zero are those defined by the equation 

G «i,'^ 2 ( x .y. u ) v ) = (V'x, Wl , Sl V'y,u )2 ,5 2 V'i W2 , S2 V'v, Wl , S i) 
By (|3.24|) . G^l (x, y, u, v) is the solution of the equation: 

i^Gsl'X 2 ) (x, y, u, v) = <5(x - v)5 W2 (y - u) - ^ 1)W2 ^ 1)S2 5(x - u)S Wl (y - v)+ 



FZ^Z\ S2 (x - y) + iC^ 4s (x - u) + Fr wi)+ (x - v)j G^ 2 (x, y, u, v) 

,(2,«) 



(3.37) 



For the two-points correlation of O x ' we are interested in the case u)\ = —uj-i = ui. For 
G1 (x, y , u, v) = G%;~$ (x, y , u, v) we find 

A+(x-y|v-y)-A+(x-u,v-u; 



G^(x,y,u,v)=e 



S w (x - v)S_ w (y - u) 



(3.38) 



Therefore, for a = C, S3 we set x = u, y = v and s = +, while for a — S1.S2 we set s = — ; 
for TC\,TC$ we set u = v, x = y and s = +; while for TG 2 , SC we set s = — . 

For the two-points correlation of Ox' a \ a ^ C, S3, we are interested in the case u)\ = w 2 = oj. 
If (x, y, u, v) = G";" a , (x, y, u, v) we find 



A (x— y I v— y ) — A s (x— u,v— u) 



G^(x,y,u,v) =e 

- AT(x-y|u-y)— AT(x— v,u— v) 

-<>s,+e L 



S w (x- v)S w (y - u) 



(3.39) 



S w (x - u)S w (y - v) 



For a = SG we set x = y, u = v and s — — ; for a — Si , S 2 we set x = u, y = v and s = — ; for 
a = S3, C we set x = u, y = v and s = +. 

Therefore, it is easy to see, by using p. 271) and (I3.34[) . that, for |x| > 1, 



(0 (2,a) ( 2 , Q)) T = 

(o£- sc] o^ sc y 



1 



cos(2ppxY 



TT 2 Z 2 C 2 (v 2 Xq + X 2 /c 2 ) x P< t (v 2 Xq + X 2 /c 2 ) 1 "-* 

1 cos(2p F a;) v 2 x\ - x 2 / c 2 



ir 2 Z 2 c 2 (v 2 Xq + x 2 /c 2 ) 2r, p (v 2 Xq + x 2 j c 2 ) 



+ (9(l/|x| 3 ), Va 

0(l/|x| 3 ) (3.40) 



1 



v 2 x\ — x 2 /c 2 



ix 2 Z 2 c 2 (v 2 x 2 +x 2 /c 2 ) 2 ^ (v 2 x 2 +x 2 /c 2 ) 2 
where m a = 1, if a = G, Si, while m a = 0, if a = SC, TCi, and 



+ 0(l/|x| 3 ), a = S u S 2 



•?7 7 -C 7 + l/2 t = (2,G),(2,S 3 ) 

Vt ~ s(7)C 7 + 1/2 t=(2,S 1 ),(2,S 2 ) 

% + C 7 + l/2 t=(2,TC 1 ),(2,TC 3 ) 

^7 + S (7)C 7 + 1/2 t= (2,SG),(2,TG 2 ) 



(3.41) 



Let us now consider the special case g a — (i.e. r\ a = Q a — 0), which we use as a effective 
model for the Hubbard model. In this case, the equations (|3.40p imply that (Oq^Ox ) decays, 
for Ixl — > 00, as |x| _2Xt , with 



2X 4 



'2 + 2r7 p -2C P * = (2,G),(2,S,) 

2 + 2?y p + 2C P * = (2, SC), (2, TG 2 

2 + 4ry p t=(l,SG) 

2 t = (l,C0,(l,Si) 



(3.42) 
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Note that 



„ = _Ullici, f „ = J^L_ (3.43) 

2 4v p+ Vp- v p+ v p - 



Let us now define if = 2X2. c — 1 an d if = 2X2, sc — 1. By using (|3.14p . we see that 



g= (l-2^) a -^ = /(l-^4)-2, p /(l + , 4 )-2, p 
Up+Up- y (1 - 1/4) + 2^ p y (1 + z/ 4 ) + 2z/ p 

~ = {l + 2v p f-vl =R _ t 

Vp+Vp- 

An p = K + K - 2 

These equations imply that all the critical indices X t and the parameter rj p can be expressed 
in terms of the single parameter K, only depending on (72/c and (74/c. In the following section 
we will show that the coupling of the model p. II) can be chosen so that its exponents coincide 
with the Hubbard ones; then, by some simple algebra, one can check the validity of the scaling 
relations (|1.22p . 

3.6 Fine tuning of the parameters of the effective model 

Let us call Vh = (g2,h,])4:.h, o"h)i h < 0, the running coupling constants in the effective model with 
ultraviolet cutoff 7^ and parameters 

91,1. = , g\\=gx = 92,n , 94 = <?4,jv , S = S N (3.45) 

so that, in particular, c = vf(1+Sn), and put vn = (g2,N ,9i,N , Sn)- We call Vh = (<?2,/i, <?4,fi, &h), 
h < 0, the analogous constants in the Hubbard model, while Vh will be defined as in £12.31 that is 
Vh = {vh,gx,hi v h)' The analysis of the RG flow given in <2]and App. |B]implies that, for h < 0, 

v h -i = v h + p(°' h) (v h , .., tJb) + r^fvH, .., v Q ,v N ) (3.46) 
v h -i = v h + ^°' h \v h , .., v ) + {v h , .., tfo, A) (3.47) 

where j3(°' h '(vh, ..,«o) is the beta function of the effective model with parameters p.45|) . modified 
so that the endpoints have scale < 0. Note that f3^°' h ^(vh, •-, i>o) is the function /3^(vh, .., vq) 
defined in (I2.4ip . modified so that, in its tree expansion, no trees containing endpoints of type 
<?i appear and the space integrals are done in terms of continuous variables, instead of lattice 
variables (the difference is given by exponentially vanishing terms). The crucial bound (IC.8[) 
and the short memory property imply that \r^ h '(vh, ■•, vn)\ < C\maxk>h |^fe|] 2 '> while the 
analysis of ^2.31 implies that r^(vh, ..,Vo,X) satisfies a bound similar to (12.511) . 

Lemma 3.1 Given the Hubbard model with coupling A such that g\$ £ D e .g, it is possible to 
choose vn as analytic function of X, so that 



92 = 2A 



0(0) - l -v{2 PF ) 



+ <3(A 3/2 ), 54 = 2Aw(0) + (9(A 2 ) , 5= 0(A) (3.48) 



and, if Vh are the r.c.c. of the effective model with parameters satisfying \3.45\ ), while Vh are the 
r.c.c. of the Hubbard model, then, V/l < 0, 

\vh-v h \<C 'g; 1 (3.49) 
1 + (a/2)|5i )0 | \h\ 

Moreover, the r.c.c. Vh have a well definite limit as N — > +00 and this limit still satisfies \3.49\j . 
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Proof - We have seen in the previous sections that the flows p. 461) and (|3.47[) have well defined 
limits v-oo and u_oo ; as /i — > — oo, if the initial values are small enough and gi.o € D e ^$. 
Moreover, the proof of this property for the flow (|3.4€>[) implies that tLoo is a smooth invertible 
function of vn, such that V-oo —vjy + 0(vfj); let us call wat(w-oo) = V-oo + 0(iP_ 0O ) its inverse. 
It is also clear that vn(v-oo) has a well defined limit as N — > oo, that we shall call v(X), and 
that this is true also for the r.c.c. Vh, h < 0. 

The previous remarks, together with f|2.73[) and (|2.74[) . imply that it is possible to choose 
vn, satisfying (|3.48[) . so that 

V^oo - V^ao = (3.50) 

In order to prove (|3.49p . we note that, because of the bound (|C.8[) and the short memory property, 
in the effective model with couplings satisfying (|3.48p . 

Iwfc-tLool <CX 2 ^ h (3.51) 

On the other hand, from Lemma 2.4 and 2.5 

h 



H-v^\<C £ [l^| 2 + A7^]<C l l q ^ (3.52) 

j= — oo 

These two bounds immediately imply (|3.49[) . | 

Let us now note that the critical indices of the effective model can be calculated in terms of 
V-oo by the same procedure used for the Hubbard model in H2.4I and that we get an equation 
like (|2.86[) . with the same function /3( '^. Hence, the above lemma allows us to conclude that the 
critical indices in the Hubbard model and in the effective model coincide, provided that the value 
of v = limN^oo ujv is chosen properly. It follows that all the indices are given by the equations 
(|3~42)1 . with 



g 2 (A) 0(0) - v(2p F )/2 3/2 

V P — "1 — A n : — = <~ U V A I 

Aire Zir sin p F 

^ 4 = MA) =A _e^ + 0(A2) 

Aire 27rsmpi? 



(3.53) 



where (13.48)1 has been used, together with c = smp F + O(X). Moreover, (]3.53f) and (|3.44l) imply 
that K = 1 - 2X[v(0) - v{2p F )/2]/(nsmp F ) + 0(A 3 / 2 ), in agreement with (TOTj) . 



4 Spin-Charge Separation 

If k ^ 0, the Fourier transform S^k + ujp F ) of the two-point Schwinger function S 2 (x) in the 
Hubbard model can be written as a tree expansion, in a way similar to eq. (2.64) of (38], whom 
we shall refer to for the notation: 

oo 

S 2 (k + p£) = £ E E E S U k ) (4-1) 

7»=0 3 'o=-ooT675 ,„, a ,o |p P^ 2 

where p F = (ujp F ,0), uj — ±. Here G 2 . ^(k) represents the contribution of a single tree r with n 
endpoints and root of scale jo; if |k| € [7 , 7 k+1 ), it obeys the bound: 

|G^(k)|<G 7 -(^o)^-E E E E ( C ' A I)" II (4.2) 

hk n=Oio=-c»T£75 ,„,k |D p £f_ v not c.p. h « 1 
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where 7^ ,n,k denotes the family of trees whose special vertices (those associated with the external 
lines) have scale or /ik+i- Moreover, d v > 0, except for the vertices belonging to the path 
connecting the root with v*, the higher vertex (of scale h*) preceding both the two special 
endpoints, where d v can be equal to 0. These vertices can be regularized by using a factor 
y-(h -jo) ^ extracted from the factor j-( h *-io)^ so that we can safely perform the sum over all 
the trees with a fixed value of h* and we get 

\S 2 (k + p F )\<C^ fl j~( hk ~ h *^ < C-t=— (4.3) 

A similar bounds can be obtained for the effective model with gx± = and couplings chosen as 
in Lemma [3.11 We shall call S£f(k) and Zh the two-point function Fourier transform and the 
renormalization constants, respectively, in this model. 
Let us put 

S 2 (k + pf) = ^-Gl(k), (k) = J-G*' M (k) (4.4) 
z '»k Z hk 

We can write 

£ 2 (k + P £) - f-SM = |Nf (k) + -^-[Gl(k) - G 2 j M (k)} (4.5) 

Note now that G^(k) differs from G 2,M (k) for three reasons: 

1) the propagators are different, which produces a difference exponentially small thanks to 
(|2.39[) . (|2.40p and the short memory property; 

2) the r.c.c. Vh and Vh are different, which produces a difference of order gx,h k i thanks to 
(I3.49|) and the short memory property; 

3) in the tree expansion of G^(k) and of the ratios Zj/Zj-x there are trees with endpoints 
of type gxi not present in the tree expansion of G^ ,M (k) and Zj/Zj—i; this fact produces 
again a difference of order gx,h k - 

These remarks, together with the fact that there is no tree with only one endpoint in the tree 
expansion, implies that 



-L[G£(k)-G^(k)] 
Zh k 



<C\Xgx, hk \^ (4.6) 



For similar reason, we have 
o 



f^ = TT #^|^ = [l + 0(A 2 )]e O(A) ^ 51 - =[l + 0(A)]L(|k|- 1 )°W (4.7) 

Zh * jit Zh ^ Zh ° 

where L(t), t > 1, is the same function defined in Theorem ll.il 

Theorem 1.3 easily follows from (I4.6[) . (|4.7p and the explicit expression (|3.34p of S u (x) in the 
effective model, applied to the case g a = 0, c = vp{l + 5). 



5 Susceptibility and Drude weight 

The effective model is not invariant under a gauge transformation with the phase depending 
both on u> and s, if gx,± > 0; however, it is still invariant under a gauge transformation with the 
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phase only depending on uj. This is true, in particular, if the interaction is spin symmetric, that 
is if g\\ = g± — 3i._i_, see item d in App. [Bj Since also the Hubbard model is spin symmetric, it is 
natural to see if one can use this "restricted " gauge invariance to get some useful information on 
the asymptotic behavior of the Hubbard model, by comparing it with the effective model with 
Si.JL > 0. 

Let us put g± = 32, 3i_l = 3i and g\\ — 32—31- We want to show that we can choose 
the parameters of the effective model gi, 32, 34, <5, so that the running coupling constants are 
asymptotically close to those of the Hubbard model. This result is stronger of the similar one 
contained in Lemma |3. 11 since now all the running couplings are involved, and this implies also 
that the values of §2, 34 and 5 are different with respect to the analogous constants defined in 
Lemma \'S. II The main consequence of these considerations is that we can use the restricted WI 
of this new effective model to get non trivial information on some Hubbard model correlation 
functions, not plagued by logarithmic corrections. 

Let lh — (gi,h,g2,h,g4:,h,Sh), h < 0, be the running coupling constants appearing in the inte- 
gration of the infrared part of the effective model. The smoothness properties of the integration 
procedure imply that, in the UV limit, lo is a smooth invertible function of the interaction param- 
eters I — (31,32,34,^), whose inverse we shall call l(lo)', hence we can fix the effective model by 
giving the value of Iq and by putting I = 1(Iq). In a similar way we call cjh = {gi,h, 32, /i! 34, ft, 5h), 
h < 0, the running couplings of the Hubbard model with coupling A. 

We now define Xh = lh — g*h, h < 0. By using the decomposition (|2.4ip for g^ and the similar 
one for we can write 

x h -i =x h + $?tf h , ..,30) - fth\lh, To)] + ffi\g h , u h ...g , u , A) + ffi\T h , .., f , f) (5.1) 

where coincides with the function f3( h ) defined in (|2.41[) . In the usual way, one can see that 



|A| + sup \l k \ 



k>h 



E^^l^l (5-2) 



k=h 



and that |/3^ 2) | < C\X\^ h , \/3 h 3) \ < C[sup fc >^ \l k \} 2 . Note that the different power in the coupling 
of these two bounds is due to the terms linear in A in the beta function for Sh, which are present 
in the Hubbard model, while similar terms are absent in the effective model, see remark after 
(1B~91) in App. [Bj 

We want to show that, given A positive and small enough, it is possible to choose Iq, hence 
xo, so that x-oo = 0; we shall do that by a simple fixed point argument. Note that a?_oo = if 
and only if 

h 

% = - E {0h 1 \9 h ,--,9o)-Pl 1 \l h ,-Jo)]+Ph 2) +P h 3) } (5-3) 

h— — 00 

We consider the Banach space Aid, 1} < 1, of sequences x = {xh}h<o with norm ||x|| = 
sup fc<0 \xk\"f~^^ 2 ^ k and the operator T : M§ — > M$, defined as the r.h.s. of (|5.3p . Given 
£ > 0, let = {x e M.{> : \\x\\ < ^A}; if A is small enough, say A < £0 an d ^A < e , and 
lh = <)h + Xh, the functions fyj^{lh, •-, h) and /3^ 3 '(i/j, .., Iq, f) are well defined and satisfy the 
bounds above, even if x is not the flow of the effective model corresponding to Iq. Hence, we 
have: 

7 -WV h \T(x) h \ <c A(£A + l) J2 7 |fc < ciA(l + £A) (5.4) 

k— — 00 
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so that i3j is invariant if £ = 2ci and A < e% — min{£ , e /(2ci), l/(2ci)}. Moreover 

h 

T(x) h - T(x') h = J2 {$h\{& + *kh>h) ~ + ? k }k>h)] 

h— — oc 

+ 0h\{9k + x' k }k> h ) - $H{9k + Zk}k>h)}} 

and \T(x)h — T(x')h\ < C2 A|]a: — x'\\, thanks to the fact that all the terms in the r.h.s. of this 
equation are of the second order in the running couplings. It follows that, if C2A < 1, T is 
a contraction in B^, so that (|5.3p has a unique solution x^ in this set; moreover, if we put 
lh = 9h + %h 1 {lh}h<o is the flow of the effective model corresponding to a value of I such that 

\9h ~ lh\ < C\X\^ h (5.5) 
Finally, this solution is such that I is equal to 170 at the first order; hence, by using (|2.26l) . we get 
g x = 2Xv(2p F ) + 0(A 2 ) , g 2 = 2Xv(0) + 0(A 2 ) , g 4 = 2Xv(0) + 0(X 2 ) , 6 = 0(A) (5.6) 



Thanks to the bound (|5.5[) . this choice of I, allows us to extend to the Hubbard model 
Lemma 1 of [27] . proved for the spinless fermion model. Hence, we can say that there are 
constants Z = 1 + 0(A 2 ), Z 3 = 1 + 0(A) and Z 3 = v F + 0(A) such that, if k < 1 and |p| < k, 

n c (p) = Z 2 3 (p p p- p ) (9) + A c + (p) 

^(p) = -z 2 (3 p 3_ p ) (9) + a j + (p) 

where (•) denotes the expectation in the effective model satisfying (|5.5f) . A c and A,- are suitable 
O(l) constants and 

Px = sV'x.w s ix = ^ ^x^ s ■ (5.8) 

Moreover, if we put p F = (cop F , 0) and we suppose that < n < |p|, |k'|, |k' — p| < 2k, < # < 1, 
then 

Gf (k' + p£, k' + P + p F ) = Z 3 (p p ^^+ +p J 9) [1 + 0(«*)] 
Gf(k' + p£,k' + p + p£) = ^(jpV^Vi+pJ^Il + 0(«*)] (5.9) 
S 2 (k' + p£) = (^^i u J i9) [l + 0{^)} . 

where G 2,1 (x) and G 2,1 (x) are defined after (|1.13[) . while the functions (/3 P V>k',u;V'k'-|-p J) and 
(ipV'k'^V'k'+p coincide with the functions (|3.16[) and (I3.I7[) . respectively, with c = v F (l+S). 
As already mentioned, if g± > 0, the effective model is still invariant under a spin-independent 
phase transformation; hence the WI (13.6[) is satisfied, if we sum both sides over s and we substitute 
Vs(p) with 

v?(p) = {5 u ,i[5 s ^ig 2 + 5 sA (g 2 - gi)] + 4,, -1^.-154} , c = v F {l + 6) (5.10) 

47TC 

Therefore we get the a WI similar to (|3.15l) , that is 

- ipo[l - P 4 (p) - 2P P (p)](/5 p ^ k ^^+ +p ^) (9) + cp[l + P 4 (p) ~ 2^(p)]( Jp ^ k ^^+ +p ^) (9) 

= I [<VWVW> (9) - ^k+P.a^k+P,^)^] (5-H) 
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where 



MP) = 9^, ^(p) = ^|^(p) (5.12) 



By replacing (15. 9[) in (|5.11l) . and comparing with (|1.14l) we get, if i^(0) = ^4, ^ p (0) = ^ 



■y = (1 ~ v * - 2^p) , — 



(1 - Di - 29 p ) , -± = c(l + P 4 - 2P P ) (5.13) 



( C) (c) ( S ) = 1 -^(l-P 4 + 2^) + y(l + P 4 -2P p ) 

1 P P/ 7rZ 2 C«2 + p 2 + c 2 v 2 p 2 V*V 

/Ac)Ac)s(g) _ 1 -pg(l-P 4 -2Pp) + c 2 p 2 (l + P 4 + 2i/p) 



^ (P) = + p4 ~ + ^ l + 2Dp)] ^TW? + 0(p) 

£(p) = + ' 4 + 2pp) + _ * 4 " 2 ' p)] ^^w + ° (p) 



(5.14) 



Moreover, by proceeding as in derivation of (|3.19[) . we get: 
Hence, by some simple algebra, we get: 

i^j" = ^/-' pll g+ ( y Up)1 + <-> < 5 - 15 ' 

where 

*V = Vp-/v p , + , = (1 - M^4) 2 - 4P 2 (5.17) 
Therefore the charge and current density correlations are given by: 



(5.18) 



From the WI (|1.14p we see that 

n c (0,p )=Q, D(p,Q) = Q (5.19) 
and this fixes the values of the constants A c and Aj in (|5.7|) . so that 



(5.20) 



If we insert (|5.13p in the previous equations, we get, for the susceptibility (11.9[) and the Drude 
weight (|1.12p . the values 

- - _ 2 _/ 1 + P 4 - 2P P + v?(l - P 4 + 2P P = 



(5.21) 

„ c(l + VA—2v n ) r/ . o, ., Kcv n 
D = -i _ 2 4 _ 2 W 1 + P4 + 2P P + ^1 _ p 4 _ 2v p )\ = p - 



nV p,+ V p 



where 





1/4) 


-2D p 




z/ 4 ) 


+ 2D p 



K= y p > -± = . ) \ ) (5-22) 
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so that 



D c 2 v p 



and this completes the proof of Theorem 11.21 

Remark We are unable to see if cv p = vp(l + S)v p coincides with the velocity cv p = Vf(1 + 
S)v p appearing in the two-point function asymptotic behavior (|3.34[) . with v p given (see f|3.14[) ) 
by 



((1 + v A ) 2 - 4v 
((1 - Av*) 



v p = ih b ( 5 - 24 ) 



with v p and defined as in (13.53[) . In fact, it is easy to see that S is equal to S at the first 
order and this is true also for v p and v p by (|5.6p . (|5.12[) and (|3.53p ; however, our arguments are 
not able to exclude that the two velocities are different. Moreover, (|5.22[) and (|3.44[) imply that 
K = K at first order, but they also could be different. Note that the equality of K and K, would 
imply that n = K/v, with v = cv p being the charge velocity, a relation proposed in |17) which, 
together with (|1.22p and (JT724J), would allow the exact determination of the exponents in terms 
of the susceptibility and the Drude weight. 

A The gi map 

Let us consider the following map on the complex plane: 

9n+i = 9n - a n gl (A.l) 

where a n is a sequence depending on go, such that, if \go\ is small enough, 

a n =a + <7 n , \<j n \ < c \g \ , (A. 2) 

for some positive constants a and Co- We want to study the trajectory of the map (|A. 1|) . under 
the condition that 

go e D StS = {z eC : \z\ < e,\Axg (z)\< it -S}, 5 €(0,77/2) (A.3) 

We shall first study the properties of a sequence g n , which turns out to be a good approximation 
of g n . Let us define: 



_^ n— 1 

A n = - V a k (A.4) 



k=0 

Lemma A.l Given 5 € (0,7r/2), there exists £q(5) such that, if e < £o(S) and go € D e _s, the 
sequence 

1 + g Q nA n 

at any step n > is well defined and does not exit the larger domain D eij s l; for £\ = 2e/(sin5) 
and S\ = 8/2. 

Proof - First of all, we choose e so that 

c e < a/2 a/2 < »a„ < 3a/2 , |3a„| < c l^o | (A.6) 

where Co is the constant defined in (| A.2|) : we can write 

A n = a n + ifi n , On > a/2 , \f3 n \ < c a \go\ ■ (A. 7) 
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Define z n := 1 + g nA n := 1 + gona n + w n ; then, if g G D e ,5, 



m i - f . c sin <5 

|1 + g na n \ > maxj smd, — — (1 + \g \na n ) 

In fact, it is trivial to show that |1 + g^na n \ > sin<5; on the other hand, if l^ol^n > 2, 

|1 +5o«a™| > |ffo|«an - 1 = (|ffo|«Q!n + 2|so|na„ - 3)/3 > (|so|raa n + l)/3 
By using (IA.8|i . we get 



(A.8 



|1 + gona n \ a sin 8 



It follows that, if e is small enough, 



\%n\ > ^ sin<5 



so that, in particular, the definition (IA.5I) is meaningful. 

Now we want to prove that g n € D eit s lt with E\ — 2e/(sin 8) and 8i — 8/2, if e 
enough. Let 50 = Poe iS °; by using (|A.8[> and (|A.9|1 . we see that, if e is small enough, 



< 



2M 



< 



2c 



\^ + a n g n\ sm.8 



(A.9) 

(A.10) 
is small 

(A.ll) 



besides it is easy to see that 

9o 



Arg 



1 + a n g n 



Arg 



Po 



a n p n 



< \e \ <n-8 



Then, since g n 



Po 



(1 + w n ), with w n of order go, for e small enough, 
I Arg (g n )\<Tr-8/2 



(A.12) 



Proposition A. 2 Given 8 £ (0, 7r/2), i/iere exists £0(8), such that, if s < £0(8) and go € D et s, 
then 

3e <5 

g n e -D<r 2 ,<5 2 , £2 = - — j , ^2 = 7 (A. 13) 

" ~ sin ' 



Moreover, if g n is defined as in &A.5\) , 



\9n ~ 9 



,, < \g n ?' 2 



(A.14) 



Proof - We shall proceed by induction on the condition (|A.14|) . which is true for n = 0. Suppose 
that it is true for n < N; then, by using (|A. 1 1|) and (IA.12[) . we see that, if e is small enough and 

n < N, 

\9n\ < 3|?„|/2 < 3e/sin(5 , |Arg (g n )\ <ir-8/4 (A.15) 
which proves (|A.13|) . Moreover, by (|A.1[) . if e is small enough, 



Note now that 

1 1 



9n+l 9n 1 — Or, 



\9N+i\<2\g N \<3\g N \ (A.16) 



= a n + a 2 n g n + A n = -J- J- + a 2 n g n + A n (A.17) 

9n+l 9n 
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where A„ is a quantity which can be bounded by ci|g„| 2 , for some constant c%. We can rewrite 
(|A.17I) in the form 



1 



1 1 



' a n9n + A„ 



9n+l 9n+l 9n 9n 

By using (|A.6p . (|A.8|) . (|A.9|) . (|A.15|1 . (|A.16[) and (|A.18|1 . we get, if e is small enough, 



\9n+i -9n+i \ = |.9jv+i| |ffjv+i| 
9 



1 



1 



gN+i 9n+i 

M 1/2 



N 



< 3\9n\ fell Y,i 6a2 \9n\ + ri \g n f] < c 2 \g N \ 3/2 (1 + N)1/2 £ T^f 



9o I 



71=0 
< 



(l + fl5o|iV) 
where C2 and C3 are two suitable constants 



(A.18) 



(A.19) 



B Symmetries of the Effective Model and RG Flow 

The RG analysis of the effective model will be done by exploiting some symmetry properties of 
a more general model, obtained by adding to the interaction (|3.3[) the term g^V^(ip), with 

The integration of the positive (ultraviolet) scales N, N — 1, .., is essentially identical to the 
one described in §2 of or §3 of [HS], for the spinless case; it does not need any localization 
procedure. 

The integration of the infrared scales is done in a way similar to the one in the Hubbard model 
described in However, before starting the multiscale IR integration, we have to perform some 
technical operations, which will make possible to compare the flow of the running couplings with 
that of the Hubbard model. 

After the integration of the UV scales up to j = 1, the free measure propagator is given by 
<7p'2(x), defined as in (|3.2p with N = 1. In this expression, the velocity c has the role of the 
Fermi velocity vf of the Hubbard model. In order to match the asymptotic behavior of the two 
models, we can not choose c = Up) for this reason we introduced the parameter 8. However, 
it is not possible to compare the RG flows of the two models if the two velocities are different; 
hence, we have to move from the free measure to the interaction the term proportional to 6. 
Moreover, since also the cutoff function x^'^iK^o, cfc)|) depends on (5, we have to "modify" it in 
X[Li](\(ko,v F k)\). 

The simplest way of performing these operations without introducing spurious singularities 
is the following one. We start with a free measure of the form 

P(d^) = AT 1 exp J £ a(k)[-iko + uv F (l + *)A]&£X£i ] } (B.2) 

where C z (k) = X [ *' 1] (IOo, ck)\)-\ We can move to the scale 1 effective interaction the term 

-*§ E "WScK = - 6VsW ' with 

k,oJ,s 
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Vs(ip) = J dx ^£,"A iuJV F d ^x^,s (B.3) 

U1,S 

The new free measure differs from (|B.2|) because S is multiplied by it;(k) = 1 — X (K&Oj cfc)|). 
On the other hand, since <5 will be chosen of order A, x''' 1 ' (\{ko, ck)\) and x''' 1 ' (|(&0j t>F&)|) differ 
only for values of k of size 7 or 7' and 

uj(k)=Q, ifx [/+1 ' Ol (l(fco^Ffc)|)>0 (B.4) 

Hence, we can write 

X [ '' 11 (|(A : o,cA : )|)=x (1) (k)+X [ ' +1 ' 0l (l(^,^fc)|)+X ( ' ) (k) (B.5) 

with x s ) an d X s ) smooth functions, whose support is on values of k of size 7 or 7', 
respectively; moreover, if we define 

Q(k)= [x [ ' +1 ' 01 (l(^o,^fc)|+X (0 (k))] _1 (B.6) 

then Ci(k) = 1, if 1 > \{k ,v F k)\ > -y l+1 and Ci(k) _1 x ( ' ) ( k ) < 1- :t follows that the free measure 
P{d,Tp^- L ^) can be written as P(d'tij^)P(d'tp^' ^), where fo 1 ' is a field whose covariances has the 
same scale properties of ip^i while 

P(d^) = N- 1 exp J - J2 ^(k)[-i*d + w«f(1 + «/(k)5)fc]^S^,S I (B.7) 

1 k,a;,s J 

The integration of the single scale field t/>D can be done without any problem. At this point, we 
start the multiscale integration of the field by performing the effective potential localization 

and the free measure renormalization as in the Hubbard model. Thanks to the support properties 
of u/(k), the steps from j = to j = l + l will give the same result we should get if the propagator 
of ■i/;!'' ] were equal to 

1 1 y^ c »kx Xi,o(IOo^Ffc)l) 
Z L 2 ^— ' — iko + ujvpk 

k 

This means that the renormalized single scale propagator will have the form (|2.38p , corresponding 
to the leading behavior of the single scale propagator in the Hubbard model. This property is 
not true only in the last step, j = I, but this is not a problem, since we have to study the RG 
flow at fixed j and I — > — 00 and, moreover, the contribution of the IR scale fluctuations to the 
Schwinger functions at fixed space-time coordinates vanishes as I — > — 00. 

Let us now analyze in more detail the RG flow of the effective model for h < 0. The main 
difference with respect to the Hubbard model is that (|2.24[) has to be replaced by 



V^\^ j ^) = g 1 , ± , j F 1 , ± (^Z' j ij) + gyF ]l (^Z-^)+ 



where the functions F a (ip) are defined as the functions V a {ip) of (|3.4[) with h(x — y) = <5(x — y); 
the absence of local terms proportional to ?p + ?p~ is a consequence of the oddness in k of the free 
propagator. The running couplings verify equations of the form 



9a,h-i ~ 9a,h = (gh,Sh, -9o, S ,g,d) 
dh-i - S h = B s h) {g h ,5h, ..ga,S( h g,d) 



(B.9) 



where a = 1 _L,||,_L,3,4 and gj = {gi,±j, g3j, g\ij, g±j, gi,j)' Note that the functions Ba 
and Bg are of the second order in their arguments; in the case of B s h \ this follows from the 
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structure of V(ip) (see p.41) 1 ). which does not allow to build Feynman graphs of the first order 
in g. For the same reason 

5 = 5+O(e 2 ), e =mzx{\g\,\5\} (B.10) 

and this relations can be inverted, if Eq is small enough. 

There are some symmetries which is important to exploit. For notational simplicity, we will 
write (G li ±,G 3 ,G\\,G±,G4, A) or (G, A) in place of (g h ,6 h , ..g ,S ,g,S) . 

a. Spin U(l). Free measure and interactions are invariant under the transformation 

where a s is a spin-dependent angle. This means that the only possible local effective 
interactions must have as many ip+ as ip~ , for each given s. Therefore, all the allowed 
local quartic interactions are the ones listed in (|3.4[) ; it is also clear from the symmetries 
u — > — lo and s — > — s that they must occur in the same linear combinations. 

b. Particle-hole symmetry. The free measure is invariant if the particle hole switching involves 

the spin s' only, i.e. ^ us -> '/Cwk w s' interactions are not invariant and we find: 

B[ h [(G h± , G 3 , G,|, G ± , Ga, A) = -Sf \~G 3 , -G 1>± ,G h -G x , -G 4 , A) (B.ll) 
-B^^Gi^x, G3, G||, Gj_, G4, A) = G3, — Gi,_l, G||, — Gj_, — G4, A) "1 

J B^ ) (G 1 ^,G 3 ,G||,G ± ,G 4 ,A) = -S^(-G 3 , -Gi,_i_, G|| , -Gj_, -G 4 , A) > (B.12) 
(Gi,±, Ga, G||, Gj_, G4, A) = — B[ h \— G3, — Gi,j_, G|| , — Gj_, — G4, A) J 

5^* (Gx,j_, G3, G||,Gj_,G4, A) = b£\— G3, — Gi,j_, G||, — G^, — G4, A) (B.f3) 

c. Chiral U(l). The free measure is invariant under the transformation 

fo,u, s -+ e isa "$L, u , s (B.14) 

for a chirality-dependent angle. All the interactions are U(l) invariant, but for V3; 
if 33 = 0, then an interaction V3 won't be generated by the flow. This fact can be seen 
also graph by graph. Indeed, in the graphs for Bu, B±, B^ and B$, the number of the 
half-lines ip+ has to equal the number of the half-lines ip~ (regardless the spin label); this 
can only happen when there is an even number of interactions V3. Therefore these three 
beta functions are even in G 3 and then, by (IB.12|) and (|B.13|) . also even in G\ t ±\ hence, if 
a = ||, JL, 4,5, 

BW(G 1>± ,G 3 , G h G ± ,G A , A) = &W(Gl ± , G 2 3 , Gp G±, G 4l A) (B.I5) 

where G^ denotes the tensor {g a ,j9a,j'}j.j'>h- By a similar argument, -B3 has to be odd 
in G3 and B^ even; then, by (|B.11I) . B[ h ^ is also odd in G\ and B^ even, so that 

B^ h ± (Gi t ±, G3, G||, Gj_, G4, A) = Gi t i_B^{G\ ± , G3, Gm, Gj_, G4, A) 

(h) - 2 ' 2 ( ' B ' 16 ' ) 

B 3 {Gi y ±, G3, G||, Gj_, G4, A) = G 3 B(h)(G 3 , G x ± , G\\ , — G±, — G 4 , A) 

where G a Ba is a shorthand for J2j>h dajB^ 1 '^. 

In this way we have found two invariant surfaces in the space of the interaction parameters 

(g,s): 

Ci = {9, S : gi t j_ =0} C 3 = {g, 5 : g 3 = 0} 
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d. Spin SU(2). It is convenient to rewrite the interaction as 

Vty) = 9i,± {Vi,xW>) - V n ty)) + (.9,1 + g h± - g ± )V\\ 

+ g x (V ± ^j) + V^j))+g 3 V 3 ^)+g 4 V 4 (^)+5 h V s (^) 

It is evident that V\ t ±_ — Vji, V± + Vji , V 3 V4 and Vs, as well as the free measure, are invariant 
under the transformation of the fields 

s' s' 

for U £ SU(2). While Vji isn't: if gu + gi — g± = it will remain zero. Thus we find four 
other invariant surfaces: 

= {9, S : gx,± = g± - g\\} , G 3 ,+ = {g, 5 : g 3 = g± + g\\} 

Ci - = {9, S : -gi,± = g± - g\\} , C 3 _ = {g, S :-g 3 = g± + g\\} 

(in fact Ci-, C 3 ,+ and C 3 ,_ are obtained form C lt+ through (|B.11[) . (|B.12I) . ()B.15|) and 
(IB. 16[) ) . C 3: + is also called Fowler invariant (see [T2], page 220). 

e. Vector-Axial Symmetry. All the terms in V(tp), but Vi{ip) and Vs^), are invariant under 

the transformation 

therefore the surface gi ,j_ = 53 = is invariant. 
Finally we consider the flow of Zh and the renormalization constant z9~' associated with 



•h 
mo 

fore, that Z h ^/Z h = 1 + b{ h \g, A) and Z^/zj? = 1 + B { p h) (G,A), with b£\g,A) 

II 



the density operator p XM . s — i/jx US i// xus in the generating functional (|3.1I) : is defined 



as z9~' G ' in (|2.27[) . It is easy to see, by using the symmetry properties of the model as be- 



B ( a h) {G 2 l7 G 2 3 ,G h G ± ,G 4 ,A) fora = z,p. Hence 



Z W 7 (i) 

- -^'1 + B^ h \G,A)} (B.18) 



with 

B (h) (G,A) = sf ) (G^ ± ,G^,G||,G ± ,G4,A) (B.19) 

C Vanishing of the Beta Function 

We recall the main ideas of the proof of (j2~4"5l) and ([2~%5Jl (see [22], [53], [23], [33] for more details). 
We consider the model (|3.1[) with = 0; 8, gu and are small but arbitrary parameters. We 
take the limit N — > 00 at fixed I; if |k| = 7' (so that, in particular, /;(k) = 1) 

S w (k) = <^>+^) ; - z,A„(k) [1 + 14/2(0(14)1 (c - 1} 

where ( • ), denotes the expectation with propagator (|3.2I) . 

A(k) = -ifc + wv F (l + Si)k 

and 

\wP(k)\ <C(s 2 +ga m ) (C.2) 
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with gi = max)i>i max{|g2.h|, |<74. h\} an d = maXj>j max{j 3 -, |<5j|}, Note that we have included 
in -Dz(k) a correction of the free Fermi velocity due to the local term 5iV$, so that w| (k) does 
not contain terms of order except those associated to irrelevant terms. Moreover, by the 
analogue of (112) of [22] 

(Pp,»Ju, u JZ- p , u , s )i = Z, (1) &(k)&(k _ p )[l + wg( k )] (C.3) 

where |W 2 ( ^(k)| < C*(ef + &7*)- In we have written the WI directly in the / — > — oo limit, 
but they are true also when I is finite, up to corrections studied in detail in §4 of [22] and §4 of 
[23]; using the bounds is such papers for the corrections we derive from the finite I analogue of 

(EH 

Dn(p)(p p ,uJ^JL p ,u, s )i = &(k-p) - &(k) + £>iv(p)flz(k ) p) (C.4) 
where k and k p are of size -f l , 

-Djv(k) = — iko + tuvp(l + S)k 
and ^ 

|i2 I (k,p)|<C 7 - 2I ^-( e ?+s«7 W ) 

Hence, if we put k = k and p = 2k, with |k| = ■y 1 , we get D N (\t)[z\ X) /Z{\ = A(k)[l + A(k)], 
with I A(k) < C(ef + ga m ), which implies that 5i = 6 + O(ef) and z[ 1] jZ x = 1 + O(ef). On the 
other hand, the value of Sh is independent of the infrared cutoff, if h > Z+l, and |<5;+i — <5;| < Cef. 
It follows that, for any j g [l,N] and uniformly in the cutoffs, 

-2^ 



*,-=* + 0(eJ) (C.5) 
(i) 

l + 0(e?) (C6) 



z 3 

In the same way, we can also use the the Schwinger-Dyson equation combined with the 
Ward Identities with a finite infrared cut-off; we get p. 241) with a correction term, which can be 
bounded as in §4 of [23]. If we restrict the resulting expression to the four point function with 
momenta at the infrared scale, we get, see [23], [21] 

9\\,j = .9|| + 0(e, 2 ) , g±,j = g x + 0{e ]) , =g 4 + 0{e 2 3 ) (C.7) 

Let us call ba Q?||> <?X) <?4> <0 the function which is obtained by B„ {gh,5h, ■ ■9o^o,g,S) (defined 
in (|B.9|) L by subtracting the contribution of the trees containing endpoints of scale grater than 

and by putting (gj,5j) — (0,0,g«,g±,gi,6), Vj = h, (the value of the r.c.c. is independent 

of the scale). Then, by Lemma 3.4 of [33, we get the bound: 

b { JHg\\,9±,9A,S)\ < C[nmx{ |.9|||, I.9 J, M,|,5|}]V J , a = \\,±,4,8 (C.8) 

In we define in a similar way the function b^(g^,g±,gi,8) in terms of B^ (gj,5j, ..go,So,g,S) 
(defined in (|B.19|> ). (IC.6|) implies the bound 

b ij \g\\,9±,94,S) <C[m a x{|.9|||,| 5± |,| 54 |,|<5|}]V J (C.9) 

The bounds (|C8[) and (|C.9|) allow to show that, in the model (|3.1j) . the infared cut-off can be 
removed by an analysis very similar to the one in Ej2.3H2.4l 
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Let us now consider the Hubbard model. In (|2.41l) we have written its Beta function as sum of 
two terms, the second of which is asymptotically negligible, by (|2.42l) ; the first term, denoted in 
(|2.41[) by pa (gj,5j] go, Sq) = /3« [G\, G 2 , G4, A) (using a notation similar to that used after 
(|B.9p ) coincides with the Beta function of the effective model on the invariant surface Ci.+ n C3, 
if we subtract from it the contribution of the trees containing endpoints of scale grater than 0. 
Hence, by using (|B.15|) and (|B.16|) . we get 

iU)m. n» n. A\ — n.aU) 



j3 1 (Gi, G 2 , G 4 , A) — Gi/3 1± (G 1 , 0, G2 — Gi, G2, G 4 , A) 
^ j) (G u G 2 , G 4 , A) = /8^(G?,0, G 2 ~G 1 ,G 2 ,G 4 , A) 
tiP(G x ,G2, G 4 , A) = pi j) (Gl 0, G 2 - d, G 2 , G 4 , A) 
fis (Gii G 2 , G4, A) = (G^, 0, G2 — Gi, G2, G4, A) 



(CIO) 



where (Gf , 0, G 2 - Gi, G 2 , G 4 , A) denotes the value of (Gf , 0, G 2 - Gi, G 2 , G 4 , A), after 
the subtraction of the trees containing endpoints of scale grater than 0. Therefore, if a ^ 1, 
the contributions of order and 1 in Gi of /3^' (G, A) are the same as the contributions of the 
same order of 0(O,O,G 2 - Gi, G 2 , G 4 , A). By using (|C.8p . one sees immediately that these 
contributions are of order 7 15 3 as functions of j, so that (|2.45[) is proved. In the same way, using 
(fCJl) and (lB~T9)) we can prove (223)- 

This completes the proof of the boundedness of the RG flow in the Hubbard model, as well 
as in the effective model with g it ± > and g\\ = g± — gi : ±. Note that, in order to prove the 
boundedness of the flow of the spin-symmetric Hubbard model, we need information from a non 
spin symmetric model; in fact, we have derived (|2.45[) from the model (|3.ip with g± ^ g^ and 
91,1. = 0. 
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